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We study the two-dimensional attractive Hubbard model using the mapping onto the half-filled 
repulsive Hubbard model in a uniform magnetic field coupled to the fermion spins. The low-energy 
effective action for charge and pairing fiuctuations is obtained in the hydrodynamic regime. We 
recover the action of a Bose superfiuid where half the fermion density is identified as the conjugate 
variable of the phase of the superconducting order parameter. By integrating out charge fiuctuations, 
we obtain a phase-only action. In the zero-temperature superconducting state, this action describes 
a collective phase mode smoothly evolving from the Anderson-Bogoliubov mode at weak coupling to 
the Bogoliubov mode of a Bose superfiuid at strong coupling. At finite temperature, the phase-only 
action can be used to extract an effective XY model and thus obtain the Berezinskii-Kosterlitz- 
Thouless (BKT) phase transition temperature. We also identify a renormalized classical regime of 
superconducting fiuctuations above the BKT phase transition, and a regime of incoherent pairs at 
higher temperature. Special care is devoted to the nearly half-filled case where the symmetry of the 
order parameter is enlarged to S0(3) due to strong q = (■7r,7r) charge fiuctuations. The low-energy 
effective action is then an S0(3) non-linear sigma model with a (symmetry breaking) magnetic field 
proportional to the doping. In the strong-coupling limit, the attractive Hubbard model can be 
mapped onto the Heisenberg model in a magnetic field, which reduces to the quantum XY model 
(except for a weak magnetic field, i.e. in the low-density limit of the attractive model). In the 
low-density limit, the Heisenberg model allows to recover the action of a Bose superfiuid, including 
the (Vp)^ term (with p the density), and in turn the Gross-Pitaevskii equation. 

PACS numbers: 71.10.Fd, 74.20.Fg, 05.30.Jp 



I. INTRODUCTION 

Many superconducting systems such as high-Tc super- 
conductors, organic conductors, heavy fermions systems, 
as well as ultracold atomic Fermi gases cannot be under- 
stood within the BCS theory. For instance, in high-Tc 
superconductors, the low dimensionality reinforces the 
role of phase fiuctuationsji*^ The short coherence length 
in these systems also suggests that they might be in 
an intermediate regime between the weak-coupling BCS 
limit of superfiuid fermions and the strong-coupling limit 
of condensed composite bosonsiii In ultracold atomic 
Fermi gases that are now experimentally available, it ap- 
pears possible to monitor the evolution from the BCS to 
the Bose limit 

A superconducting system at low temperature is conve- 
niently described by a low-energy effective action written 
in terms of a few relevant (bosonic) variables. The mini- 
mum description requires to consider the phase of the su- 
perconducting order parameter, but other variables such 
as the amplitude of the order parameter or the charge 
density may also be included. The effective action (or 
the corresponding equations of motion) is sufficient to 
describe macroscopic quantum phenomena such as the 
Meissner effect, the fiux quantization, the Josephson ef- 
fect, or the vortex dynamicsi^ In a Bose superfiuid, the 
low-energy effective action leads to the Gross-Pitaevskii 
equation^Sii i.e. a non-linear Schrodinger equation for the 
complex (superfiuid) order parameter * where — p 
is the condensate density. In a Fermi system, there 



is in general no simple relation between the amplitude 
of the superconducting order parameter and the den- 
sity. In the strong-coupling limit, where fermions form 
tightly bound pairs which behave as bosons, we expect 
the Gross-Pitaevskii equation to hold. Moreover, since 
Fermi and Bose superfiuids should behave similarly in 
many respects, a Fermi superfiuid should be described by 
a non-linear Schrodinger equation similar to the Gross- 
Pitaevskii equation even in the weak-coupling Hmiti^ 
This is the conclusion reached by previous works j^°i^^i^^ 
although the "wavefunction" in this case is not the su- 
perfiuid order parameter (except in the strong-coupling 
limit). 

The low-energy effective action provides a convenient 
framework to discuss the BCS-Bose crossover between 
the weak-coupling BCS limit and the Bose limit of pre- 
formed pairs In two dimensions (2D), 
there is an additional motivation to introduce an effec- 
tive action written in terms of the phase of the order 
parameter. The Berezinskii-Kosterlitz-Thouless (BKT) 
phase transitionii2*2fi which takes place in a 2D Fermi 
superfiuid is clearly out of reach of fermionic approxima- 
tions based on diagram resummations like the T-matrix 
approximationii 

While the 2D attractive Hubbard model is not an ap- 
propriate microscopic model for most superconducting 
systems of interest, it can be used to understand a num- 
ber of general issues relevant to many cases. The main 
characteristics of this model are well-knowni2i*2^ Away 
from half-filling, there is a BKT phase transition to a 
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low-temperature superconducting phase. Long-range or- 
der sets in at zero temperature and breaks SO (2) symme- 
try. At half-filling, q = pairing and q — (tt, tt) charge 
fluctuations combine to form an order parameter with 
SO (3) symmetry. The superconducting transition then 
occurs at zero temperature and breaks SO (3) symmetry. 
In the weak-coupling limit, superconductivity is due to 
strongly overlapping Cooper pairs, and Bogoliubov quasi- 
particle excitations dominate the low-energy physics. In 
the strong-coupling limit, fermions form tightly bound 
pairs that Bose condense at low temperature thus giv- 
ing rise to superfluidity. Pair-breaking excitations are 
not possible at low energy and the thermodynamics is 
controlled by (collective) phase fluctuations. 

In this paper, we derive the low-energy effective action 
for charge and pairing fluctuations in the 2D attractive 
Hubbard model on a square lattice. We discuss the BCS- 
Bose crossover and the phase diagram. We use the map- 
ping of the attractive model onto the repulsive half-filled 
model in a magnetic fleld which couples to the fermion 
spins. In this mapping, the charge and pairing flelds 
transform into the three components of the spin-density 
fleld. At half-flUing, the magnetic fleld vanishes in the 
repulsive model. This case has been studied in detail 
in Refs. 12^241 At zero temperature, the system evolves 
from a Slater to a Mott-Heisenberg antiferromagnet as 
the interaction strength (i.e. the on-site Coulomb repul- 
sion) increases. Because of the SO (3) symmetry of the 
order parameter, the Neel temperature vanishes in agree- 
ment with the Mermin-Wagnes2^ theorem. At flnite tem- 
perature, the system is in a renormahzed classical (RC) 
regime with an exponentially large antiferromagnetic cor- 
relation length. Away from half-fllling, the magnetic fleld 
reduces the symmetry to SO (2) so that a BKT phase 
transition occurs at flnite temperature. In the attractive 
model, this corresponds to the suppression of q = (tt, tt) 
charge fluctuations at low energy. The main advantage 
of studying the repulsive model is that the SO (3) sym- 
metry of the order parameter at or near half-fllling can 
be easily handled. 

In Sec. ^ we map the attractive model onto the re- 
pulsive one by a canonical particle-hole transformationSfi 
and then obtain the effective action for the spin fluc- 
tuations in the presence of a flnite magnetic fleld (i.e. 
away from half-flUing in the attractive model). Col- 
lective bosonic flelds are introduced by means of a 
Hubbard-Stratonovich decoupling of the Hubbard inter- 
action. This crucial step in our approach differs from 
the usual decoupHng in two respects. First, we write the 
interaction in an explicit SO (3) spin-rotation form by in- 
troducing a unit vector fir at each site and timeiSSiSiSLSS 
This allows to recover the Hartree-Fock (HF) theory at 
the saddle-point level, while maintaining SO (3) spin- 
rotation symmetry (in the absence of the magnetic fleld) . 
Second, we introduce two auxiliary real flelds. One is 
simply the Hubbard-Stratonovich fleld m^^ for the am- 
plitude of the spin density at site r. The other one, m^, 
is directly connected to the actual amplitude of the spin 



densitypS^ While these two flelds are proportional at the 
saddle-point (i.e. Hartree-Fock) level, they differ when 
fluctuations are taken into account. In the weak- and 
strong-coupling limits, ^Jlr can be identifled with the 
spin density, which allows a direct interpretation of the 
low-energy action S[m, fl] in terms of physical quantities. 
Because of the magnetic fleld, the spin component along 
the fleld, S"^, takes a flnite value and its fluctuations are 
small at low energy. Since spin amplitude fluctuations are 
also small, the important fluctuations correspond to rota- 
tions around the magnetic fleld axis. We derive the effec- 
tive action S[m, fl] in this case. S[m, fl] takes a simple 
form in the weak-coupling (Slater) and strong-coupling 
(Mott-Heisenberg) limits. These two Hmits differ in the 
role of the Berry phase term. 

In Sec. Ill CA we deduce the effective action S[p, A] of 
the charge (p) and pairing (A — |A|e'®) fluctuations in 
the attractive model. By integrating out amplitude fluc- 
tuations (|A|), we recover the action S[p,Q] of a Bose 
superfluid where half the fermion density is identifled as 
the conjugate variable of the phase 6 of the supercon- 
ducting order parameter. This action is parametrized 
by the mass rrih of the "bosons" and the amplitude g of 
the repulsive interaction between "bosons", mt and g 
are computed as a function of particle density and in- 
teraction strength. We then analyze in more detail the 
weak (BCS) and strong (Bose) coupling limits. In the 
BCS limit, we flnd that the action is a function of charge 
and phase fluctuations only, since amplitude fluctuations 
of the superconducting order parameter decouple. In the 
strong-coupling limit, the amplitude of the superconduct- 
ing order parameter and the fe rmion density p satisfy the 
relation |Ar| = ^^/pr{2 — p^). The Bose superfluid ac- 
tion ^[p, 6] = 5[|A|,9] then entirely describes the dy- 
namics of the superconducting order parameter A. In 
Sec. Ill Dl we derive the phase-only action ^[O] by inte- 
grating out charge fluctuations. S[Q\ corresponds to an 
0(2) sigma model with an additional term proportional 
to the flrst-order time derivative of 8. The phase stiffness 
and the velocity of the 0(2) sigma model are obtained 
as a function of particle density and interaction strength. 
At zero temperature, superconducting long-range order 
gives rise to a gapless (Goldstone) phase mode smoothly 
evolving from the Anderson-Bogoliubov mod^ at weak 
coupHng to the Bogoliubov mode^^ of a Bose superfluid 
at strong coupling. In Sec. Ill El we show how we can ex- 
tract from the phase-only action an effective (classical) 
XY model whose phase stiffness is a function of density, 
interaction strength and temperature. This allows us to 
determine the value of the BKT phase transition temper- 
ature as a function of density and interaction strength. 
The XY model also yields an estimate of the crossover 
temperature Tx below which the system enters a RC 
regime of phase fluctuations. At higher temperature, for 
Tx < T < Tpair, there is a regime of incoherent pairs 
[Cooper (local) pairs at weak (strong) coupling] with no 
superconducting short-range order. Tpair is estimated 
from the HF transition temperature. 
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In the vicinity of half-filling, the analysis of Sec. ^ is 
not sufficient since q = (tt, tt) charge fiuctuations (in the 
attractive model) are not considered. For a weak mag- 
netic field (in the repulsive model) , there are strong fiuc- 
tuations of S"^, and the analysis of Sec. |n] breaks down. 
This case is dealt with in Sec. IIIII We show that the 
dynamics of spin fiuctuations in the repulsive model is 
governed by an SO (3) non-linear-sigma model with a 
(symmetry-breaking) magnetic field proportional to the 
doping. The magnetic field defines a characteristic tem- 
perature Tso(3^2) above which the S0(3) spin-rotation 
symmetry is restored. Below Tso(3^2)) the system en- 
ters a RC regime of spin fiuctuations with SO (2) sym- 
metry. The global phase diagram, as a function of den- 
sity, interaction strength and temperature, is discussed 
in Sec. Ill F| [For clarity, we discuss the phase diagram at 
the end of Sec. and postpone the technical analysis of 
the SO(3)^SO(2) crossover near half-filling to Sec.lml] 

In Sec. IIVI we consider the strong-coupling limit in 
more detail. First we show that the repulsive Hubbard 
model reduces to the Heisenberg model to leading order 
in \/U. This allows to obtain the collective modes be- 
yond the long-wavelength approximation. We then show 
that except for a strong magnetic field (i.e. in the low 
density limit of the attractive model), the Heisenberg 
model in a magnetic field reduces to the quantum XY 
model. In the low-density limit, we recover from the 
Heisenberg model the usual action of a Bose superfiuid 
(for bosons of mass 1/J— U/At'^ and density p/2), in- 
cluding the terms proportional to (Vpr)^ that were omit- 
ted in Sec. ^ The classical equation of motion derived 
from this action is nothing but the Gross-Pitaevskii equa- 
tion. We therefore obtain a correspondence between the 
Gross-Pitaevskii equation in the attractive model and the 
semiclassical spin dynamics in the repulsive model. 

To our knowledge, there is no systematic study of the 
2D repulsive Hubbard model in a magnetic field. When 
translated in the language of the attractive model, our 
results reproduce, in a unique framework, a number of 
previously known results. We also obtain new results, 
in particular regarding the BCS-Bose crossover and the 
phase diagram. 



II. AWAY FROM HALF-FILLING 

The attractive Hubbard model on a square lattice is 
defined by the Hamiltonian 

= -^4(£+^)Cr - t/^nrt^r^ (1) 



where t is the nearest-neighbor hopping operator: 

— i(Cr-f-x ~t- Cr— X "l" Cr-\-y ~^ ^r — y)- (2) 

X and y denote unit vectors along the x and y axis. With 
the notation of Eq. 1^, —U is the on-site interaction with 



L'^ > in the attractive case. The operator cj^ (Cro-) cre- 
ates (annihilates) a fermion of spin a at the lattice site 
r, Cr — (cr|,Crx)"^, and rira — cl^Cra- At half-filHng, 
particle-hole symmetry implies that the chemical poten- 
tial fi equals —U/2. In the following, we will consider 
only hole doping so that /i < —U/2. We take the lattice 
spacing equal to unity and h = = I. 

Under the canonical particle-hole transformation^^ 



Crl - (-1)^ 



(3) 



the Hamiltonian becomes (omitting a constant term) 



t + /i + 



u 

2" 



U 



r 

(4) 

where (cr"^, cr^, cr^) denotes the Pauli matrices. The trans- 
formed Hamiltonian Q corresponds to the repulsive half- 
filled Hubbard model in a magnetic field 



ho = fi + 



U 
2" 



(5) 



along the z axis coupled to the fermion spins. In the at- 
tractive model, the chemical potential ^ is fixed by the 
condition (cjcr) = po = 1 — x where po is the mean 
density and x the doping. Under the particle-hole trans- 
formation 



c\a^'Cr -I- 1. In the repulsive model. 



the magnetic field /lo is then determined by 

(cJcr^Cr) = -X. 



(6) 



The charge-density and pairing operators transform as 



Pr-4cr ^ {2Si + l) 



(7) 



where = cl^Cr (i' = x,y,z) and = ± iS^. 
In the repulsive model, spin fiuctuations are clearly the 
collective (bosonic) fiuctuations of interest. Accordingly, 
in the attractive model, both pairing and charge fiuctu- 
ations should be considered on equal footing. One of the 
motivations to study the repulsive model is that one has 
to consider only the particle-hole channel. On the con- 
trary, a direct study of the attractive model would require 
to consider both the particle-hole and particle-particle 
channels in order to take into account the charge and 
pairing fiuctuations. The simultaneous introduction of 
auxiliary Hubbard-Stratonovich fields in these two chan- 
nels is not without problem if one requires the saddle- 
point approximation to recover the HF (or mean-field) 
theoryji^ Furthermore, in the attractive model at half- 
fiUing, q = pairing and q = (tt, tt) charge fiuctuations 
combine to form an order parameter with SO (3) sym- 
metry. The SO (3) symmetry of the order parameter is 
much more easily handled in the repulsive model, where 
the distance from half-filling determines the (symmetry- 
breaking) magnetic field /iq (see Sec. IIII|I . 
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We can write the partition function of the repulsive 
model as a path integral over Grassmann fields c*^,Cra, 
with the action 



(8) 



where t is an imaginary time and /3 — 1/T the in- 
verse temperature. i?[c^, c] is obtained from the Hamilto- 
nian Q by replacing the operators by the corresponding 
Grassmann fields. We now introduce auxihary bosonic 
fields for the collective spin fluctuations in a way that 
fulfills the three following requirements: i) the HF (or 
mean-field) approximation is recovered from a saddle- 
point approximation; ii) the SO (2) (or SO (3) if ho = 0) 
spin-rotation symmetry is maintained; iii) the auxihary 
fields correspond to the spin-density field Sr ~ cl^Cr not 
only at the saddle-point level, but also when fiuctuations 
are taken into account [cr = [u^ .a"^ 




We start from the identit\i22i2i21 



FIG. 1: Correspondence between the spin density Sr in the 
repulsive model, and the charge density Pr and superconduct- 
ing field Ar = |Ar|e*®"' in the attractive model [see Eqs. Q]. 



(4cr)'-(c 



J7r 



(9) 



In the following, we shall consider the charge field Acr 
only at the saddle-point level, i.e. Acr — i{U/2){clcr) — 
iU/2. The term —iclAcrCr in II12II cancels the chemical 



where fir is an arbitrary site- and time-dependent unit 
vector. Hr is defined by its polar and azimuthal angles 
9r,^r- Spin- rotation invariance is made explicit by per- 
forming an angular integration over fir at each site and 
time (with a measure normalized to unity). The charge 
term (cj!cr)^ is decoupled by means of an auxiliary (real) 
field Acr. In order to decouple the spin term (cJcr-firCr)^, 
we introduce in the path integral the unit factor 

1 = / V[m] Y[ S{fnr - c\a- ■ firCr) 
r,T 

- y P[TO,mHS]e-/o'''i-i:.-r(™.-4-i^.^.),(10) 

where m^^ is a Lagrange multiplier field which im- 
poses the constraint rrir = 2Sr • rir- Both m and 
TO^^ are real fields. Note that integrating out the ra^- 

field [see Eqs^^ below] one obtains the action p^o^^jnation with = mo, - „.o , ^..^ ^ . 

S{c\ c, Ac, m ^] where A and m ^ are the Hubbard- configuration of the unit vector field given by 



potential term — (J7/2)cJcr in 

Eqs. II11I12|I are the starting point of our analysis. In 
Sec. Ill Al we show that the HF theory is recovered from a 
saddle-point approximation over the auxiliary fields mr, 
TO^^ and rJr. In the following sections, we go beyond 
the HF theory and derive a low-energy effective action 
^[m, $7] for spin fluctuations, and deduce the effective 
action S\p, A] of charge and pairing fluctuations in the 
attractive model. 



A. Hartree-Fock theory 

In the presence of the uniform magnetic fleld along the 
z axis, we expect the ground state to exhibit AF order 
in the (x, y) plane and a ferromagnetic order along the 
z axis. We therefore consider a static saddle-point ap- 



proximation with rrir 



HS 



and a classical 



Stratonovich flelds which decouple the interaction term 
(El. In general, -^fir cannot be directly identifled with 
the spin density Sr, but the identification turns out to 
be correct in the hydrodynamic regime both at weak 
([/ < M) and strong {U > 4t) coupling. 

Using Q and we write the action as S* = 5'o + 'S'int 
with 



:1 



(-1)' 



sm C^qx + cos t^oz. 



(13) 



The HF action then reads 

V 



So = 

Sint — 



j^\TY,4(dr-i-hoa''^^^Cr, (11) 

rf^ r a2 



^HF = m{-^rnl 

+ / rfTVc]:[a^-f-/ia^-AHS(_i)V^]cr, 
Jo 

(14) 



-4(iAcr + im^^a- ■ J7r)Cr 



(12) 



where N is the total number of lattice sites, and 

A^^ = im^^ sin 6*0, 

h = ho + im^^ cos 00- 



(15) 
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The HF action is quadratic and can be easily diagonal- 
ized. The single-particle Green's functions are given by 

Gaik^iuj) = -(ccr(k, zu;)c* (k, iw)) 
_ -iuj - ekg 



AHS 



(16) 



where 



(17) 



Ek — — 2t(cos fci:+cos ky) is the energy of the free fermions 
on the square lattice, a = —a, Q = (tt, tt), and 
uj = TrT{2n + 1) (n integer) is a fermionic Matsubara 
frequency. C(j(k, zw) is the Fourier transformed field of 

The saddle-point equations are obtained from 
dZuF/dmo = dZup/dA^^ = dZwe/dh = 0, where Zhf 
is the partition function in the HF approximation: 



mo 
Ao 



AT 

U 
mo . 
— sm Oq 

(-l)""/ t . N 



TOO cos 6*0 = —{h-ho) 



(18) 



(19) 



(20) 



Eq. then implies 



Too cos 6*0 = Jji^~ ^o) 



(21) 



Ao = {—iy{S^) is the AF order parameter in the re- 
pulsive model, and the superconducting order parame- 
ter in the attractive model, h is an effective magnetic 
field which takes into account the mean ferromagnetic 
magnetization (cJcr^Cr) along the z axis. In the attrac- 
tive model, h — ho — xU /2 — + pqU /2 corresponds to 
the chemical potential renormalized by the Hartree self- 
energy. 

Using Eqs. ifTfiji . we rewrite the saddle-point equations 

irraimi as 



2 _ /■ tanh(/?£:kT /2) 
tanh(/3£;kT /2) 



CkT- 



■kT 



(22) 
(23) 
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FIG. 2: Chemical potential fi — ho — U/2, order parameter 
Ao ^ and 6o in the T = HF state vs doping x = 1—po for U = 
2t, 4t and 12t (solid, dashed, and dotted lines, respectively). 
The energies are measured in units of t. 



where = ^. At T = 0, Eqs. de- 

termine the superconducting order parameter A^^ (or, 
equivalently, Aq) and the renormalized chemical poten- 
tial h. Using II15I21II . we then obtain ho, im^^ and 
0Q as a function of U and x [Figs. EE). Eqs. (j22l2;i|l 
also determine the HF transition temperature T^^ where 
AF long-range order (i.e. superconducting order in the 
attractive model) sets in. At the transition, we have 
Ao = A^s ^ o+, 6*0 = TT and im}^^ = xU/2. 

Partial analytical results can be obtained in the limits 
of weak and strong couplings. 



1. Weak coupling (U <^ At) 

At half-fiUing (a; = 0), /i = /lo = and Oq = 7r/2. 
The zero-temperature order parameter and the transition 
temperature are given by 



aHS 
^0 

rriHF 



32ie 



te 



t/U 

TJu 



(24) 



6 



12 



16 




FIG. 3: Chemical potential jj. = ho — U/2, order parameter 
A^^ and Oq in the T = HF state vs U fox x ^ 0.1, 0.5 and 
0.9 (solid, dashed, and dotted lines, respectively). Here and 
in the following figures, we use the analytical results in the 
weak-coupling regime {U < t) where the numerics becomes 
difficult because of the exponentially small value of Aq^. 



Note that at half-filling, the SO (3) symmetry is restored 
and the AF order parameter can have a component along 
the z axis. The latter corresponds to a q = Q charge- 
density-wave order in the attractive model. The choice 
9q = n/2 corresponds to a state with only superconduct- 
ing order. 

Away from half-filling {Hq < 0), Bq 
xU/2 since the superconducting order parameter is ex- 
ponentially small. At T = 0, the superconducting order 
parameter is determined by 



TT and im^^ 



2 

U 



de 



A/'o(e) 

I -At 



(25) 



4t 



where E = [(e — /i)^ + 
-2/16^2)1/2]^ with ee 



e"/ Wf)^'"], witn e t [— 4t,4t], is the density of states of 
free fermions on a square lattice [K is the complete eUip- 
tic integral of the first kind]. Since the integral in Ij25|l 



is peaked around e = h for A^^ ^ 0, we have replaced 
the density of states A/q (e) by its value at the renormal- 
ized chemical potential h. In the weak-coupling limit, 
we can also neglect the effect of the order parameter on 
the chemical potential (i.e. set Aq^ = in ^^). We 
then have h ~ tp where ef is the Fermi energy of the 
non-interacting system. With these approximations, we 
obtain (for U ^ Q and x fixed) 



aHS 
^0 



2[(4t)2-e|,]i/2e-i7A/HF77^ 



_ 4t-|eF| 
2irt 

U 



if 1 - X < 1, 



(26) 



Since h ~ ep belongs to the non- interacting band 
[— 4i, 4t], the excitation gap (i.e. the minimum energy 
required to break a pair) equals 2Aq^. 

Similar arguments show that the transition tempera- 
ture is given by 



(27) 



where 7 ~ 1.78 is the exponential of the Euler constant. 
We recover the results of the BCS theory with [(4t)2 — 
e|,]^/2 playing the role of the cutoff energy [Eqs. (I26l27|l ]. 



2. Strong coupling (U ^ At) 

To leading order in 1/U (and x fixed), at T = we 
find jmg^ = U/2 (or mo = 1) and cos^o = ~x = ho/2J. 
This gives 



U , 



A- = -(l-x2)l/2, 



u 



2Jx. 



(28) 



In the strong-coupling limit, the excitation gap equals 
2(/i2 + A^s^)i/2 = f/ to leading order in 1/U. For a; < 1, 



B. Effective action S[m,ii] 

In 2D, the HF theory breaks down at finite tempera- 
ture since it predicts AF long-range order below T^^ . 
Nevertheless, the HF transition temperature bears a 
physical meaning as a crossover temperature below which 
the amplitude Aq of the AF order parameter takes a fi- 
nite value. This can be interpreted as the appearance of 
local moments perpendicular to the magnetic field and 
with an amplitude Aq — /U. Note that at weak- 
coupHng these "local" moments can be defined only at 
length scales of order ^0 ~ V^o^ ^ which corre- 
sponds to the size of bound particle-hole pairs in the HF 
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state. Thus, stricto sensu, local moments form only in 
the strong-coupling Hmit when '^-^ 1- 



Below , the fluctuations of the flelds 



,HS 



and 

around their HF values are therefore expected to be 
small. Below a crossover temperature Tx < , AF 
short-range order sets in. The AF correlation length be- 
comes much larger than the lattice spacing, and the ef- 
fective action for AF fluctuations (which correspond to 
rotations of fir about the z axis) can be derived within 
a gradient expansion. In this section, we derive the low- 
energy effective action S'[m, fl] for temperatures below 
the crossover temperature Tx- For T <^ Tx < , 
the coefflcients of the effective action (which are related 
to HF quantities) can be evaluated in the zero tempera- 
ture Hmit. This means that we neglect the exponentially 
small number of thermally excited quasi-particles which 
give rise to non-analytic contributions (Landau damp- 
ing terms) to the effective actioni22i2i As shown below, 
Tx ~ J(l — x'^)/2 ^ T^^ in the strong-coupling limit 
(J = At'^/U), but Tx ~ T^^ in the weak-coupling limit. 
Fluctuations can be parametrized by 



Srrir = 
5mf = 

Pr = 



Too, 



-,HS 



HS 



ipr - Q - r. 



(29) 



In the HF state, Srrir — Sm^° = = = 0. The 
effective action S[p,q,m^^,fn] is obtained by integrat- 
ing out the fermions, and assuming pr,5mnSm^^ and 
dfiPr,dfj,5mr,dfiSm^^,dfiqr to be small (/.t = 0,x,y and 
do = dr). We do not assume q^ to be small so that our 
approach is valid even in the absence of AF long-range 
order. It is convenient to introduce a new fermionic vari- 
able 0r = i'Pri, <t>riY' deflued by Cr = Rr(t>r where i?r is a 
time- and site-dependent SU(2)/U(1) matrix satisfying 

Rr(T ■ VLfRl = (T • rJr. (30) 

The above deflnition means that TZr, the SO (3) element 
associated to i?r, maps $7^' onto J7r- The U(l) gauge free- 
dom is due to rotations around $7^', which do not change 
the physical state of the system. The pseudofermion 0r 
has its spin quantized along TZr^. The action H11I12|I can 
then be expressed as 

S = S'hf + Si + S2 



(31) 



where S'hf is the HF action (I14II and 

51 = rfr V'^t^or^r, 

^0 

52 ^ dr ^ (0|;ylr,r'0r' +C.C.). (32) 

•^0 (r,r') 



(r, r') denotes nearest neighbors and we have introduced 



iy—x,y,z 



Ar,r' = RIRv' - 1 



E 



(33) 



We use the notation i?r = drRr- cr" is the 2x2 unit 
matrix. 

The effective action S[p, q, to] is obtained by inte- 
grating out the fermion fleld (pr and the Hubbard 
Stratonovich fleld to!?^. The flnal result reads 



S[p,q,m] = ^^[p^qJIpp{q)pq + q^qllqq{q) 



+Sm^qIl„im{q)Smq + 2p^qWpq{q)qq 
+ 2p^qlipm{q)5mq + 2q_qllqjn{q)5mq\ 

+5Sb, (34) 

where q = (q, iwi/), with cj^ (v integer) a bosonic Mat- 
subara frequency. The coeflicients 11 are given by IIB24|l 
and SSb by IIB21|I . Details of this rather long calculation 
are given in Appendix El to which we refer readers inter- 
ested in results for the repulsive Hubbard model [Eq. Q]. 
In the bulk of the manuscript, we shall focus on the at- 
tractive Hubbard model as deflned in Q. 



C. Effective action ^[p, A] 



Results of Sec. Ill Bl can be easily translated to the at- 
tractive model. Below Tpair = , the flnite amplitude 
Ao of the superconducting order parameter can be inter- 
preted as the appearance of incoherent pairs. At weak 
coupling, the size ~ ^/^o^ of these (Cooper) pairs is 
much larger than the lattice spacing (^o ^ !)• At strong 
coupling, the preformed pairs are local and are expected 
to behave as hard-core bosons (the hard-core constraint 
comes from the Pauli principle which prevents double 
occupancy of a lattice site). The crossover temperature 
Tx marks the onset of strong superconducting fluctua- 
tions (i.e. C ^ 1 with ^ the superconducting correlation 
length) . 

As already pointed out, in general ^Jlr cannot be 
directly identifled with the spin density Sr = cl^Cr- In 
order to flnd the relation between Sr and ^ftr, we add 
to the action the source term 



S.J 



■Jo 



(35) 



with Jr = (Jp , J^,Jr)- The charge and pairing flelds in 
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the attractive model are then obtained from [see Eq. 0] and H^ffl|l . we obtain 



SS 
5S 

M 

SS 
5JS 



= 2(-l)''|Ar|c0Ser, 



,7=0 



= ~2(-l)''|Ar|siner 



J=0 



J=0 



Pv - 1- 



(36) 



We then proceed as in the preceding section (see 
also Appendix El- We integrate out the fermions 
to obtain the effective action to second order in 
Pr,Smr,dfj,Pr,df_i6mr,df_tqr and first order in Jr. Using 
then Eqs. we obtain the relation between mpjJlr 

and Ar,pr: see Eq. IID6|l in Appendix IdI Eq. (ID6II takes 
a simple form in the BCS and Bose limits (see Sees. Ill (TTl 

and mm . 

To proceed further in the general case, we note that 
there is no coupling between Vg and p, q, m in the ac- 
tion S[p,q,m\ (AppendixEl- Since V0r — — V^r, and 
Q,Sp,6\A\ are functions ofp,q,m (Appendix^, S[p, A] 
takes the general form 



S[p,A] = - dr d\ 



iprQr 



(-K) 



-(VOr 



~n,p(<5pr)'+n|A||A|(<5|Ar|)^ 
-2Up\^\SprS\Ar 



(37) 



Here we use the fact that half the fermion density is the 
conjugate variable of the phase of the superconducting 
order parameter, as required by gauge invariancei^ The 
action S[p, A] was previously obtained in Ref.H^ For our 
purpose, it is not necessary to determine the expression of 
the coefficients IIpp, ni^n^i and Hpi^i which can be found 
in Ref. H^i^ To make contact with the usual description 
of a Bose superfluid, we integrate out the amplitude field 
lAI. We thus obtain the action 



Pbr = y, 



4/00 



9 = 4 Hpp- 



n 



|A||A| 



(40) 



In order to calculate g without computing Hpp, n|A||At 
and HpiAj, we integrate out charge ffuctuations. This 
yields the phase-only action 



S[Q] = dT j 



-po J d,T I d'^rQr 



(41) 



S[Q] corresponds to an 0(2) sigma model with an addi- 
tional term proportional to the first-order time derivative 
of 8. TOfc and g are related to the phase stiffness and 
the velocity c of the 0(2) sigma model: 



nib 



Pa 

„2 



(42) 



S[Q] can be directly obtained from the action S[p, q, m] 
(Sec. lII D|| . which allows to determine p^ and c, and there- 
fore nib and g. Anticipating on the results of Sec. Ill PI 
we show in Figs. I4I5I mh and g as a function of doping 
and interaction strength. Because both p° and c vanish 
in the low-density limit, the numerical determination of 
nib and g from lli2ll is difficult when x ^ 1. We have 
therefore only considered x < 0.9. 

The low-energy effective action II38I39|I ensures that 
the Fermi superfluid will behave similarly to a Bose su- 
perfluid. Indeed, from the classical equation of motion, 
we obtain the two basic equations (in imaginary time) of 
a superfluidt^ 



1 /-^ 

S[p,e] = - dr 



iprQr 



n2 



-K) 



n 



|A||A| 



(38) 



This action is similar to the superfluid action^^ 



Sb^ J dr I d'r 



ipbrQr 



Pba 
2mb 



(39) 

for bosons of mass and density Pbr {Pbo is the mean 
density and dpbr = Pbr — Pbo )- 9 is the amplitude of the re- 
pulsive interaction between bosons. Comparing Eqs. Ij38|l 



ipbr 



V • 



PbO, 



VOr = 0, 



nib 

iQr + gSpbr ~ 0. 



(43) 



For a Bose superfluid, the coefficient in front of (V0)^ 
and V9 in Eqs. ll39)l and l|i3ll is the superfluid density 
(divided by the boson mass). Here, because the coeffi- 
cients of the effective action S[p, 8] are calculated from 
the HF action, we expect the superfluid density to be 
given by the full density pbo = Po/2. In Sec. lIVI we derive 
the effective action S[p, 8] in the strong-coupling low- 
density limit [including the terms proportional to (Vpr)^ 
that are omitted in l(39jl ]. without relying on the HF the- 
ory. We obtain a boson mass nib = l/J in agreement 
with lUnil when C/ > 4i (see Eq. JHOJ below for x ~ 1). 
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FIG. 4: Boson mass mi, and interaction strength (; vs doping 
X = 1 — po for U — 2t, 4t and 12t (solid, dashed, and dotted 
lines, respectively). 
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FIG. 5: Boson mass mt and interaction strength g vs U for 
a: = 0.1, 0.5 and 0.9 (solid, dashed, and dotted lines, respec- 
tively) . 



section, we have obtained an effective action for the lat- 
tice case which is vahd for all values of the interaction 
strength. The validity of our approach is quite obvious 
both at weak and strong couplings. In the latter limit, it 
simply follows from the mapping between the attractive 
model and the Heisenberg model under the particle-hole 
transformation (pj (see Sec. lIV|l . We therefore expect our 
approach to provide a good description of the BCS-Bose 
crossover as the interaction strength increases. 

We discuss below in more detail the BCS and Rose 
limits. 



1. BCS Umit 

Using the results of Appendix IB II we deduce^ from 



5pr = —drrir 

d\Ar\ = —prX. 



(44) 



Eqs. 



2 



lUijl can be directly obtained by assuming Sr = 
• and considering the limit Aq 0. This shows 
that the identification between Sr and ^ftr holds in the 
weak-coupling limit. When inverting the particle-hole 
transformation (EJ , the kinetic energy remains unchanged 
and the HF spin susceptibility IIqq{0,0) becomes the HF 
compressibility k of the attractive model. The effective 
action S[q, m] [Eq. IIB28|l ] then becomes 



s[p,e] 



dr / d^r 



iprQr 



U 



(45) 



Amplitude fluctuations (|A|) do not couple to density 
and phase fluctuations in the weak-coupling limit. We 
also verify that half the fermion density is the conjugate 
variable of the phase 6 of the pairing fleld. The gauge 
choice = fr/mo [Eq. IIB21|l ] is crucial to obtain this 
result. For instance, with ipr = 0, one would obtain 
{pr — l)/2 as the conjugate variable of 9r. 

Prom lUfill . we conclude that in the BCS limit the sys- 
tem behaves as a Bose superfluid with pbr = Pr/2 and 



Po 



-KY 



u 



(46) 



There is nevertheless an important difference with the 
Bose superfluid case. S[p, Q] has been obtained by in- 
tegrating out amplitude fluctuations (|A|) and therefore 
does not describe the full dynamics of the superconduct- 
ing order parameter Ar = |Ar|e'®"' (except in the strong 
coupling limit, see Sees. Ill (T2l and ITV)l . 

Effective actions similar to l|38|l have been derived pre- 
viously for continuum or lattice models pi*iS*iiii^ In this 



S. Bose limit 

In the strong-couphng hmit, the auxiliary fleld 
can be identifled with the spin density Sr (see Appendix 
iDlll so that 

Pr-l = 

Ar = ^^^r", (47) 
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where fl^ = il'^iifl'^. The condition fJ^ = 1 implies that 
charge and amplitude fluctuations are not independent 
but tied by the relation 



1 



|Ar| = -v/pr(2-pr). 



(48) 



Ar = |A. 

incide 



There is therefore a one-t o-one correspondence between 
the bosonic fleld ^'r — y/pr/2e'^^'- and the pairing fleld 
In the low-density limit, both fields co- 
. since |Ar| ~ 

Prom (|17jl . we deduce Spr — — 2sin0oPr and 6\Ar\ — 
cos^oPr (with cos^o = —x). These equations can also 
be obtained from ijDfip . The effective action S[p, Q] then 
reads [see Eq. IIB? 



s[p, e] 



dr / d^r 



tPrQr+j{l~X^){VQr? 



'2J{Spr 



(49) 



Again we verify that half the fermion density is the con- 
jugate variable of the phase 6 of the pairing field. In the 
Bose limit, we therefore have 



Jil + x)' 
9 = 8J. 



(50) 



The ECS and Bose limits differ in the role of the Berry 
phase term 



(51) 



where we have included the gauge-dependent contribu- 
tion ljB20p . For U 3> 4i, the Berry phase term becomes 
|pr0r and therefore determines the entire dynamics of 
the phase 6. This is expected since the dynamics of the 
Heisenberg model entirely comes from the Berry phase 
term.^* For U ^ 4i, the Berry phase term only gives 
|po0r and does not contribute to the T = phase col- 
lective mode. The missing term, |^Pr0r, comes from the 
second-order cumulant {{Si-\-S2)'^)c- For all values of the 
interaction strength, the contribution |po dr J2r ®r to 
the action, which is responsible for the Magnus force act- 
ing a vortexjiS is given by the Berry phase term of the 
effective action S[m,n]m^ 



rP 0.1 



0.2 0.4 0.6 0. 




FIG. 6: Phase stiffness pi and velocity c vs doping x = 1 ~ 
po for U = 2t, At and 12t (solid, dashed, and dotted lines, 
respectively). 



where we have taken the continuum limit in real space. 
Eq. 1I52I1 is vaHd in the hydrodynamic regime defined by 
the momentum-space cutoff A '-^ min(l, 2Ag ^/c). Its 
validity also requires the fluctuations of 6 to be small, 
a condition which is not fulfllled in the vicinity of half- 
fllling. This case is discussed in detail in Sec. IIIII The 
(bare) phase stiffness pi and the velocity c of the 0(2) 
sigma model in Ij52|l are determined by 



-K) 



8 ' 
4 9^0 ojI 



n'„,Hs(g) 
n' Hs^Hs (9) 



n,p(g) - 



5 (-J^n^ 



n' 



3(9) 



npp(g) 



,HS('3) 



n' 



(53) 



D. Phase-only action ^[B] 

Integrating out p, 6m^^ and 5m in the action 
S'Ip, g, m^^, to], and using Qj. — —Qr, we obtain the 
phase-only action 



I dT I d\ 



(VGr)^ 



e 



"oPo / dr d rOr, 



(52) 



At zero temperature, there is superconducting long-range 
order (in the attractive model). From ll52)l . we deduce 
the existence of a gapless (Goldstone) mode with disper- 
sion u! = c|q|. This collective mode smoothly evolves 
from the Anderson-BogoHubov mode^ at weak coupHng 
to the Bogoliubov mod#l of a Bose superfluid at strong 
coupling (see Sees. Ill D IIII D 2ll . 

Figs. I6I7I show c and p^ vs U and x. Our results 
reproduce the collective mode velocity obtained from 
a random-phase-approximation (RPA) calculation about 
the zero-temperature HF superconducting state4S»^ii^ 
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FIG. 7: Phase stiffness and velocity c us [/ for a; = 0.1, 0.5 
and 0.9 (solid, dashed, and dotted lines, respectively). 



BCS limit 



In the weak-coupling limit, using the results of Ap- 



pendix |BT1 and Sec. Ill (" Tl we obtain 



U 



(54) 



We further simplify this result by noting that A[/^^/[(e- 

yields [see Eqs. llA7p in Appendix IA] 

K = 2JVo{eF) 
Using also (Appendix El 



eF)^ + AHs2]3/2 in thejimit A^s o. This 

(55) 



where Vk = Vkeic, we finally obtain 



(56) 



/kt'k'^(ek-£F) 



(57) 



is the mean square velocity on the (non-interacting) 
Fermi surface. Eqs. lf57jl reproduce the expression of the 
velocity of the Anderson-Bogoliubov mode^fi in a 2D con- 
tinuum model if we identify with vj^ = (kF/niY (with 
kp the Fermi momentum and m the fermion mass). They 
were also obtained in Ref. from an RPA calculation 
about the zero-temperature HF superconducting state. 



2. Base limit 

In the strong-coupling limit, using the results of Ap- 
pendix ^3 we obtain 

J. 



Ps 



-(1 



(58) 



We recover the velocity c — y/pbo9/n^b of the Bogoliubov 
mode in a Bose superfluid(2& where the mass mi, and the 
interaction amplitude g of the bosons are defined in Ij5()|l 
{Pbo = Po/2). At half-fiUing, x = 0, Eqs. lf58)l agree with 
results obtained directly from the Heisenberg model. 



E. Berezinskii-Kosterlitz-Thouless transition 

The effective action S[Q] derived in Sec. lIIDl is valid for 
temperatures below the crossover temperature Tx, i.e. 
T <^Tx (see Sec. lIIB|l . Since Tbkt ~ Tx ~ Tpair in the 
weak-coupling limit (as will be shown below), this effec- 
tive action is not sufficient if one is interested in the BKT 
phase transition. In this section, we derive the phase-only 
action S[Q] for temperatures T < Tpair and obtain the 
BKT phase transition temperature Tbkt- We consider 
the static (i.e. classical) limit, where non-analytic con- 
tributions due to the Landau damping teruisM^ are not 
present. In this limit, the q field decouple from other fiuc- 
tuations (p, 6m^^ and Sm), and we obtain (see Appendix 



El 



s[e] = 



Ps = 



§ I d^r(Ver)^ 

i-K) n-.(o,o) 



£k£kT 
8i?kT 



tanh ■ 



kT 



sin^ krr 



2T 



2rcosh' 



(59) 



(.60) 



The action II59|I is valid in the hydrodynamic regime and 
must therefore be supplemented with a momentum-space 
cutoff A - min(l,2A[/S(T)/c). n^^(0,0) is a current- 
current correlation function in the attractive model. It 
vanishes at zero temperature, and the phase stiffness re- 
duces to the mean kinetic energy as obtained in Sec. Ill PI 
The expression l|671|l of the finite-tem perature phase stiff- 
ness has also been obtained in Refs. l43l44l45L 

In order to extract the BKT phase transition temper- 
ature, we assume that the system is described by an XY 
model whose continuum limit is given by Ij59ll ^iiiiSiiSiiS 
Although this assumption can be justified at strong- 
coupling (except in the low-density limit) (see Sec. IIV|l . 
it is in general not correct iiS Nevertheless, it should give 
a reasonable estimate of the BKT phase transition tem- 
perature. This leads us to the 2D XY Hamiltonian 



HxY 



-P, 



{r,r'> 



COs(0r — 6r'), 



(61) 
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defined on a square lattice with spacing A~^. The XY 
model lj6T|l is known to have a BKT phase transition at 
the temperature Tbkt defined by 

Tbkt = Qp°{Tbkt), (62) 

where Q ~ 0.898 is estimated from Monte-Carlo 
simulations In Il62|l . we have written explicitly the 

temperature dependence of p^. Tbkt is obtained by solv- 
ing simultaneously Eqs. II22I23I62II . 

From the XY model, we can extract another charac- 
teristic temperature, 

Tx = 2p2(rx), (63) 

defined as the mean-field transition temperature of Hxy- 
Tx marks the onset of a RC regime of superconduct- 
ing fluctuations. In this regime, superconducting fluc- 
tuations become quasi-static, the amplitude |Ar| of the 
order parameter takes a flnite value, and strong phase 
fluctuations develop (i.e. ^ 3> 1, with ^ the correlation 
length) !^ 

In the strong coupling limit, Tbkt — Qj{l — x^) and 
Tx ~ -^(1 — x^) <C Tpair- In the weak-coupling limit, 
Tbkt - Tx - Tpair oc g-i/^^oC^F) exponentially 
small. 



F. Phase diagram 

Near half-filling, q = Q charge fiuctuations become 
important and are suppressed only at very small tem- 
perature. [This corresponds to a weak magnetic field Hq 
in the repulsive model.] We call Tso(3-^2) the crossover 
temperature above which charge fiuctuations restore the 
SO (3) symmetry of the order parameter. Postponing the 
determination of Tso(3_»2) to Sec. IIIII we discuss in this 
section the phase diagram of the 2D attractive Hubbard 
model. 

Let us first discuss the phase diagram as a function of 
doping (Figs.|HE|)- Below Tpair (not shown in the figure 
for U — 12t), the amplitude of the superconducting order 
parameter takes a finite value, but with no short-range 
order of the phase (i.e. ^ '--^ 1 with ^ the superconducting 
correlation length). This corresponds to the appearance 
of incoherent pairs (Cooper pairs at weak-coupling and 
local pairs at strong coupling) . Below Tx , superconduct- 
ing correlations start to grow. The system enters a RC 
regime of fiuctuations with SO (2) symmetry: the phase 
correlation length ^ increases and rapidly becomes much 
larger than the lattice spacing ^ 1). At Tbkt, a BKT 
phase transition occurs and the system becomes super- 
conducting. Superconducting long-range order sets in at 
T = 0. The situation is slightly different near half-fiUing 
where Tgo(3^2) < Tx- Below Tx, both superconduct- 
ing and q = Q charge fiuctuations start to grow, and 
the order parameter has an effective SO (3) symmetry. 
For T < Tso(3^2), charge fiuctuations are suppressed 
and only superconducting fiuctuations continue to grow. 
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FIG. 8: Tbkt, Tso{3-.2), Tx and Tpair vs doping a; for C/ = 
12t, 4t and It. For IJ — 12t, Tpair, which is nearly a vertical 
line around s = 1 on the scale of the figure, is not shown. 

This corresponds to a regime of fiuctuations with SO (2) 
symmetry which eventually drive a BKT phase transi- 
tion at the temperature Tbkt- Once the SO (2) regime 
is reached, the BKT phase transition occurs rapidly, so 
that Tbkt can be estimated by Tso(3-»2) (see Sec. lIIIB| . 
Therefore, the true transition temperature is estimated 
as min(TBKT, Tgo(3-,2)); where Tbkt is the transition 
temperature obtained in Sec. Ill El i.e. by neglecting the 
q = Q charge fiuctuations. The resulting phase dia- 
gram is shown in Fig. El In the intermediate and strong- 
coupling regimes (i.e. for U > At), we clearly identify a 
RC regime above the BKT superconducting phase and 
an incoherent-pair regime at higher temperatures. We 
believe the persistence of the incoherent-pair regime at 
weak-coupling (see Fig. El for U = 2t) to be likely an ar- 
tifact of our simple estimate of Tpair- Due to many-body 
effects not taken into account in our approach, Tpair is re- 
duced with respect to the HF transition temperaturei^ 
We expect Tpair — Tx at weak coupling and thus the 
disappearance of the incoherent-pair regime. We also ob- 
serve that the RC regime with SO (3) fiuctuations extends 
to higher doping at strong coupling. Fig. ^1 shows the 
phase diagram as a function of U for x = 0.1, 0.5 and 
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FIG. 9: Phase diagram of the 2D attractive Hubbard model vs 
doping X for U = 12t, 4t and 2t. SC: superconducting phase 
(algebraic order). RC-S0(2): renormalized classical regime 
with fluctuations of S0(2) symmetry. RC-S0(3): RC regime 
where the fluctuations exhibit an effective S0(3) symmetry 
due to the presence of strong q = Q charge fluctuations. FL: 
Fermi liquid. The incoherent-pair regime (IP) corresponds to 
the formation of Cooper (local) pairs at weak (strong) cou- 
pling, without short-range SC order. Superconducting long- 
range order sets in at T = 0. 



0.9. 

Away from half-filling, the BKT phase transition tem- 
perature derived from the criterion Tbkt p?(Trkt), 
with p° given by ijfiCijl . was also discussed in Refs. 143. 
|4'4.45..4fL Good quantitative agreement with Quantum 
Monte-Carlo simulations was obtained4^ Our estimation 
of the crossover temperature Tx for U = At also a gree s 
with the Quantum Monte-Carlo simulations of Ref. 

In the weak-coupling limit, we expect a pseudogap 
to appear in the density of states and the spectral 
function ^(k, = — 7r~^ImG(k, w) below the crossover 
temperature Tx [G(k, w) is the single-particle Green's 
function] This pseudogap results from strong scat- 
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FIG. 10: Phase diagram of the 2D attractive Hubbard model 
vs U for x = 0.1, 0.5 and 0.9. 



tering of particles on collective fluctuations in the RC 
regimei^iiS In the BCS limit, the pseudogap region is 
small, except in the vicinity of half-filling where an SO (3) 
RC regime is observed over a wide temperature range 
(0 < T < Tx) [see Fig.Eli^ As the interaction strength 
increases, the weak-coupling pseudogap should progres- 
sively evolves into a strong-coupling gap due to the for- 
mation of local singlet pairs when T < Tpair- Pseudogaps 
in the 2D attractive Hubbard model have been seen both 
in Quantum Monte Carlo simulations^ii^Si^ and analyt- 
ical approaches At half-filling, the transfor- 
mation of the weak-coupling pseudoga p into the strong- 
coupHng gap has been studied in Refs. l23l24 



III. NEAR HALF-FILLING 

In this section, we calculate the crossover temperature 
Tso(3^2) introduced in Sec. Ill FI The approach of SecHTI 
holds only when q = Q charge fluctuations (in the attrac- 
tive model) are weak for all temperatures below Tx- In 
the repulsive model, this corresponds to the limit where 
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the magnetic field hg is strong enough to suppress fluc- 
tuations of the 5^ component of the spin density. Near 
half-filling, q = Q charge fiuctuations are suppressed 
only at very low temperature. We therefore expect a RC 
regime with SO (3) symmetry in the temperature range 
Tgo(3-,2) <T <Tx followed by a RC regime with S0(2) 
symmetry at lower temperature (Tbkt <T < Tso(3^2))- 

In the repulsive model, a weak magnetic field ho = 
/iqz can be treated perturbatively. When it vanishes, 
the low-energy effective action for spin fluctuations is a 
NLcrMiS2i2iiS& From general arguments?^ we expect the 
magnetic field to modify the kinetic term of the NLcrM, 
(Ar — 2iho X rir)^, as obtained in the Heisenberg 
limit 1^ [Here n (n^ = 1) is the Neel field describing 
low-energy AF fiuctuations.] We show below that this 
is indeed correct. We consider the low-temperature limit 
T ^ Tx where the coefficients of the effective action can 
be evaluated in the zero-temperature limit. 

Following Ref.E3. we introduce only one auxiliary field 
(mr) to decouple the spin term (cJ(T-J7rCr)^ in ©i^The 
action then reads 



order, we obtain 

S = S'hf + Sp Si + Sd + Si2 + Shg + Ssn 
5rrii + 2moSmr 



Sho 



U 



Sri 



Sa = ^ 







dr 



E 



s, = 



Sl2 — 



-mo 



dr 



C.C. 



'ho 



dT^(l)lRla''Rr 



dr 



t — hoa^ 



rf-k \ ' r 

mr<T ■ IZrjCr + -jj- 



(64) 

At the HF level, neglecting the magnetic field, we have 
TTir = Too and Clf = (— l)'"z. We choose the AF mag- 
netization along the z axis as in Ref. il The saddle- 
point equation reads 2mo/U — (— l)'"(cJ(T^Cr) . Follow- 
ing Haldanei^fiifii in the presence of AF short-range order 
(T < Tx), we write 



Or = (-l)'"nrVl -L?+Lr 



(65) 



rir is the slowly varying Neel field, whereas Lr is a cant- 
ing vector, orthogonal to rir, taking account of local fer- 
romagnetic fiuctuations. We assume Lr to be small, 
which clearly restricts the validity of this approach to 
weak magnetic fields. As in SeclUl we introduce a new 
field (f) defined by Cr = Rr(t>r where Rr is a time- and 
site-dependent SU(2)/U(1) matrix satisfying 



RrCr^'Rl = cr • rir. 



(66) 



TZr, the SO (3) element associated to Rr, maps z onto 
rir. We also define the rotated canting field lr = 7?.~^Lr. 
Given that 7?."^ rir = z and Lr -L rir, the lr vector lies in 
the {x — y) plane. 

In order to express the action in terms of the field, 
it is convenient to make use of the SU(2) gauge field 
A„ 



T,^=x,y,z Kr^'' defined as 

^Or = —RldrRr, 



(67) 



Since the gauge field is of order 0{d^), we can expand 
the action with respect to 1, A^, hg and 5m. To second 



Ssm — 



(68) 



where the spin-density current j^^ is defined in (IA2|I . For 
ho — 0, S reduces to the action derived in Ref. 24. S'hf is 
the HF action. Sp and Sd are paramagnetic and diamag- 
netic terms, respectively. Si and S12 are first-order and 
second-order corrections in 1. Ssm is the contribution due 
to amplitude fiuctuations.^^ 

The effective action ^[n, L, Sm] is obtained by integrat- 
ing out the fermions. To second order in A^, 1, 6m and 
ho, one finds 

S[n,L,Sm] = S[n,L; ho = Sm 0] + {Sp + Sho + Ssm) 



+ 2{Sho + Ssm){Sp + Sl))c 

\ ^ Smi ~\- 2mo5mr 
dT^- 



U 



(69) 



where the averages (• • • ) are taken with the HF action, 
^[n, L; ho = 5m ~ 0] is the action with no magnetic 
field (and no amplitude fiuctuations) and was derived in 
Ref. I24I {Sp) is a Berry phase term. It was ignored in 
Ref.|2J, since it does not play any role in the RC regime of 
the /lo = NLcrM. In order to calculate the HF averages, 
we write 



Sho 



'ho 



/3 



dr 



(70) 



where Bj^, is defined by Rla'R^ = j:.=.,y,. B^r^'' . A';, 
and Bor are calculated using 



Rr = 



(I) 
in(%) 



cos 
sin 



62 



(71) 
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Here we define rir by its polar and azimuthal angles 9^ 
and ifr- The angle ipr comes from the U(l) gauge freedom 
in the definition of Rr- We obtain 



where 



-K) 



AX 


^ ^^rSinV^r- 

2 


I 

-ify smt 
2 


?r COS Ipy, 




i ■ 

= -Oy COS ?/;r + 


i 

-ify sin 


9y sin ?/;r , 


4^ 

^Or 


i 

= —Ipr cos Oy + 








— — sin Oy cos i/i 


ri 






— smOySinipy, 








— cos 9y. 







(72) 

^Jlr (m 7^ 0) is obtained from Aq^. with the replacement 
dr -idfj.. 

We find (see Appendix El 



(Sp) 

(Sho) 



-i-ff I dT^{-lY{ipyC0Sey+1^y), 







0, 



{Si 



{ShoSi) 



Jo 

dr^ho • (rir X hy), 
mongj / dr Vho • Lr, 

"'0 



(73) 



and (ShoSsm) = (SsmSp) = (SsmSi) = 0. There is there- 
fore no coupling between amplitude and direction fluc- 
tuations in the limit of a small magnetic fleld. The HF 
correlation function ngg = Hpn (0, 0) is given in Appendix 
lA 21 Using the result of Ref. ,24, for S[n, L; ho = Sm ^ 0], 
one has 



/9 

S[n,L] = I dr 



i-K) 
16 



(Vnr)^ 



H — 22. (hr - 2iho x rir 



yvxx 



--mongoLr • (rir X hr - 2zho) 
+5*5 [n], 



(74) 



where we have taken the continuum limit in real space. 
(K) is the mean kinetic energy in the HF state (for 
ho — 0). We denote the Berry phase term {Sp) by S'_B[n]. 
Integrating out the L fleld with the constraint Lr -L rir, 
we flnally obtain 



S[n] = ^ 



dr / d 



(Vnr) 



2 (Ar — 2iho X rir)^"' 



(75) 



-K) ( 1 



TJXX 



(76) 



PI and c are the spin stiffness and the spin-wave velocity 
in the absence of fleld (/iq — 0, i.e. at half-fllling in the 
attractive model) . It can be checked analytically that the 
expression of c agrees with lf53jl evaluated at /iq = 0. The 
action iTZfill is valid in the hydrodynamic regime deflned 
by the momentum-space cutoff A min(l, 2mo/c). In 
the strong-coupling limit, c = v^J and = J/4, we 
recover the NLtrM obtained from the Heisenberg model 
in a magnetic fleldi^ The crossover temperature Tx can 
be obtained from the criterion ^ ~ 1, where ^ is the AF 
correlation length deduced from the NL(jM»S4 However, 
since Tx is weakly doping dependent near half-fllling, 
we can also consider the estimate obtained in Sec. Ill Fl 
[Eq. (jSHl]. 



A. Low-temperature limit 

In the low-temperature limit, one expects to recover 
the results obtained in Sec. |n| Let us flrst consider the 
NLcrM Il75p within a static saddle-point approximation 
where the Neel fleld rir' li^s in the {x,y) plane. The 
classical action reads 



mi 



(77) 



The magnetic fleld is determined by the condition 
{Npy^dha In^' = -X. From (jnl, we deduce 



[cla^Cy 



ho = —X 



(78) 



The chemical potential ^ = ho — U/2 obtained from l|78|) 
is in very good agreement with the result of Sec. ^ for 
X < 0.2 (Fig. EJ- In the strong-coupling limit, c = \/2J 
and p° = J/4, we flnd ho = —2Jx as in Sec. Hll 

Let us now consider fluctuations about the Neel state 
deflned by rir'iS^ We have shown in Sec. Ill (^1 that the 
phase of the superconducting order parameter is ob- 
tained from Or = -arctan(f7^/17J) -f Q • r [Eq. ^ ]. 
From l(65|l . we then deduce Qy — —ify. It is also clear 
that the 9 fleld is related to the q = Q charge fluc- 
tuations. In the low-temperature limit \n^\ <C 1. The 
low-energy effective action S[Q,9] is obtained by assum- 
ing IVGrl, lOrl, |V6lr|, |^r|, \9y " 7r/2| < 1. For the cal- 
culation of the Berry phase term S'B[n], we choose V'r = 
{—iy(pyU /(2too). As in Sec^ this gauge choice ensures 
that electromagnetic gauge invariance is satisfledi^ This 
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FIG. 11: Zero-temperature chemical potential /i = ho — U /2 
vs doping X for U = 2t, At and 12t. The dashed lines show 
the result obtained from Eq. Il78|l . The solid lines reproduce 
the results obtained by solving Eqs. I|22l23ll [see Fig. |2]. 



^ives S[Q,0]^ Sci + S[Q] + S[0] where 



e 



s[e] 



2 7o 



-2P0 / dT I d?Q, 



I dT I d\ 



2" 



(79) 



(80) 



5 [9] corresponds to the ho = limit of the phase-only 
action Ij52|l derived in Sec. |n| S[6] describes q = Q 
charge fluctuations with a gap 2\ho\. These fluctuations 
were not considered in Sec. ^ We show in the next 
section that they play an important role when \ho\ is 
small, i.e. near half-fllling in the attractive model. 



B. High-temperature limit: SO(3)^SO(2) 
crossover and BKT transition 



The preceding results are valid at low temperatures 
when fluctuations of 9^ are small. In this regime, the 
SO (3) rotation symmetry is broken by the magnetic fleld. 
At higher temperatures, the gapped mode will be ther- 
mally excited thus restoring the SO (3) spin- rotation sym- 
metry. In this section, we determine the crossover tem- 
perature Tso(3^2) separating these two regimes and the 
BKT transition temperature Tbkt < 7so(3^2)- 

The SO(3)^SO(2) crossover can be understood from 
renormalization group (RG) arguments. Following the 
standard procedure,'^, we write the Neel field rir as ((1 — 
nr^)i/2 Hr), Hr = (n^,n^), with iHrl < 1. This yiclds 
the spin-wave action 



dT 



d^r 



TT2 4/;2 



Si, 
(81) 



with momenta bounded by the NLcrM cutoff A. In agree- 
ment with Sec. nil A1 we find a gap in the (bare) propaga- 
tor of the field. The interaction part Si of the action is 
evaluated aX — Q. RG equations are obtained by inte- 
grating out degrees of freedom with momenta between A 
and Ae"''' and rescaling momenta, energies and fieldsiS^ 
At the beginning of the RG procedure, we can ignore the 
gap in the fluctuations of 11^ and treat the magnetic fleld 
perturbatively. We then obtain the usual RG equations 
of the SO (3) NLcrM (/iq = 0) for the dimensionless cou- 
pling constants g = cA/pP and t = T/p°, together with 
the flow equation dho/dl — hg. The RG flow should be 
stopped when we reach the strong-coupling regime of the 
S0(3) NLcrM {t{l) - 1) or when the magnetic field can- 
not be treated perturbatively anymore {\2ho{l)\ ~ cA). 
Introducing the characteristic length A~^e', the former 
condition defines the SO (3) AF correlation length ^ while 
the latter defines the magnetic length £,ho = c/|2/io|. 
We are therefore left to consider two different cases. If 
C ^ ^ho, t{0 ~ 1 occurs before |2/io(0l ~ cA. The system 
is then disordered by SO (3) fiuctuations before reaching 
the SO (2) regime. This will occur above the crossover 
temperature Tso(3-2) defined by C(7so(3^2)) ~ £.ho-— H 
C/io ^ i-e- T < T'so(3-2), an SO(3)^SO(2) crossover 
takes place. For |2/io(OI > cA, fiuctuations of 11^ are sup- 
pressed, and only fiuctuations of 11^ survive. This regime 
is described by an SO (2) NLcrM with coupling constants 
t{lho) and g{lho) where Iho is defined by |2/io(Z/iJ| cA. 
The very existence of the S0(2) regime impHes t{lhf,) < 1. 
The latter inequality is also approximately the condition 
for the SO (2) NLcrM to be in the low-temperature BKT 
phase. This implies that the temperature range where 
the system is disordered by SO (2) fiuctuations is very 
narrow, so that we can identify the BKT transition tem- 
perature with the SO(3)^SO(2) crossover temperature: 
Tbkt ^ Tso(3— 2)- 

In order to determine rso(3-.2), we use the following 
expression for the ho = correlation length ^i^i 



C(T)=A^^exp(^ 



(82) 



where K ~ 0.05 and ps is the zero-temperature spin stiff- 
ness in the Neel state. To estimate ps, we use the one- 
loop RG result^ 



Ps 



1 



cA 
47rp0 



Tso(3^2) is therefore given by 

27rps 



Tc 



SO(3-*2) 



2Trps 



In 



\2K\ho\ J \ \ho\ J 



(83) 



(84) 



where we have taken K = 0.05. 730(3^2)) obtained from 
l(84ll . is shown in Fig. |H| 

The definition of 7so(3^2) is meaningful only below 
the crossover temperature Tx which marks the onset of 
AF short-range order and defines the temperature range 
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where the NLcrM holds. Once 730(3^2) becomes of the 
order of Tx, as \ho\ increases, q = Q charge fluctuations 
(in the attractive model) can be ignored. In this regime, 
the analysis of Sec. Ull holds. 

IV. STRONG-COUPLING LIMIT 



where we have used {RlRr')aa{RlRr)aa = (1 - ' 

ftr')/2. We have introduced the spin— i coherent 
state irir) = Rr\ T) cos(|^)e-5(vr+^r)| |) + 

sin (^) e^^'^'"'''''^! ipr is arbitrary and corresponds 

to the U(l) gauge freedom in the definition of the SU(2) 
matrix Rr- We therefore recover the action 



In the strong-coupling limit U 3> we can directly 
integrate out the fermions to recover the action S[^^\ of 
the Heisenberg model in a mag netic field (Sec. WJl . 
We can then go beyond the hydrodynamic limit (q 
0) considered in Sec. O and obtain the collective mode 
dispersion over the entire Brillouin zone. When |/io| is 
weak enough (i.e. away from the low-density limit in the 
attractive model), the Heisenberg model reduces to the 
quantum XY model fSec. lIVBl . In the low-density limit, 
the Heisenberg model allows to recover the usual action 
of a Bose superfiuid, including the terms proportional 
to (Vpr)^ that were omitted in Sec. ^ and in turn the 
Gross-Pitaevskii equation fSec. lIVC)l . 



A. Heisenberg model 



S= ["^ drj^ [(Orljir) - ho • fir] + J XI 

(89) 

of the Heisenberg model in a magnetic field ho . 

Consider first the classical ground-state defined by 
rj^' = (— 1)'' sin 6'ox + cos 6'o2;. Minimizing the classical 
action 



-ho cos Oq + — 003(26*0) 



(90) 



with respect to 60, we find 



We write the action Ijllll2ll in terms of the (j) field 
defined by Cr = Rr4>r and Rrcr^Rl — cr ■ J7r: 

f/9 



^0 
COS0O = ^ 



(91) 



S = S^t+ [ dT\^(l)l{RlRr - hoRla'' Rr)(f>r if l^ol < 2J and 



t {(t)lRlRr" 



+ C.C.) 



(85) 
(86) 



where we neglect spin amplitude fiuctuations and ztoq ^ = 
U/2 to leading order in 1/U (see Secnj. S^t is the action 
in the atomic limit {t = 0). The effective action of the 
angular variable fl is obtained by integrating out the 
fermions. To lowest order in t/U and hg, we obtain 



TT otherwise (for ho < 0). The 
ci.a'-Cr) — —X translates into (fi^) = —x (see 
Appendix ID 111 , i.e. cos 6*0 = —x within the classical ap- 
proximation. We thus obtain ho = —2Jx as in Sec. lII A 2l 

The Heisenberg model allows to obtain the collec- 
tive excitations of the attractive Hubbard model in the 
strong-coupling limit without taking the continuum limit. 
We introduce the variables pr = {Or — 0o)/2 and qr — 
(fir — Q ■ r and derive the effective action to quadratic 
order in Pr and qr — qr' (for r, r' first neighbors): 



S = 



S = I dTY,[{<f>lRlRr<l>r) ' hoic^lRla'RrCl^r)] 
r 

^ "'0 ri,r2 



: sin 6'o(Pr9r - Pr^r) + 2/10 COS BqpI 



'^Paqr 



sin^ 60 



+ 4 E 

{r,r') 

{qr -qr'Y 



-2cOs{2Bo){pr+Pr'f 



(87) 



(92) 



where try equals t for r,r' nearest neighbors and van- 
ishes otherwise. The averages in l(87jl are taken with the 
atomic action. They can be easily calculated using the 
parameterization Ij71|l of i?r. One finds 

{c^lRlRr(j)r) = (f^rlrir), 

{(jjlRlo'' Rr(t)r) = ^r, 

{ {(t>l Rl Rr, C^rM^' RI' Rr', ^r^r' ) 



5r{,rJr',,r,5{T - r')jj{l - ^r^ ' ^r,), 



The Berry phase term (rJrl^^r) has been evaluated with 
ipr = Vr- This gauge choice is similar to the one made in 
Sec. ni We thus obtain 



S 



'-P0 f dTY,qr + \Y.^P-^^^-i)'^''^i) ( ' 
J, g \ / 

(93) 
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FIG. 12: Collective mode dispersion ujq in the strong-coupling limit as obtained from the Heisenberg model [Eq. 19511 ] for 
X = 0.1 (solid line), 0.5 (dashed line) and 0.9 (dotted line). [F = (0,0), M = {tv,tt) and X = (7r,0).] 



— 4Jcos(26'o)(l + 7q) + 4/io cos 6*0 i^i, sihOq 

— w,ySin6'o Jsin^ f?o(l — 7q) 



(94) 



where 7q — {cosq^ + cosqy)/2. Eqs. II92I94|I assume the fluctuations of p to be small and are therefore valid only for 
ho ^ 0, i.e. away from half-filling in the attractive model. At zero-temperature, there is AF long-range order (i.e. 
superconducting order in the attractive model). Collective modes are then obtained from detl?~^(i3') = with the 
analytic continuation to real frequencies iw^ — > Wq. This gives^ 



= 2hlil~^^)-i2hl-AJ^){l-^l) 

= Sj'x'il - 7q) - 4j2(2x2 - 1)(1 - 7^). 



(95) 



Fig- El shows ujq for different values of the doping x. For q — > 0, we obtain a spin-wave mode (satisfying p ~ 0) with 
dispersion ujq ~ c\q\, c — \/2J\^l — x^. This mode corresponds to the Bogoliubov mode obtained in Sec. ITnT2l In 
the vicinity of Q = (vr,7r), we find a mode with the dispersion 



(2/io)'±V2J|l-3x2|i/V 



(96) 



where the + (— ) sign refers to the case x < 1/VS {x > l/\/3). When x < l/\/3, this mode corresponds to a local 
minimum of the energy with a gap |2/io|- It involves fluctuations along the magnetic fleld axis which correspond in 
the attractive model to q = Q charge-density fluctuations. At the critical value of the doping Xc = l/\/3, this local 
minimum becomes a local maximum. Note that the value \2ho\ of the gap was also found in Sec. IIII Al for a weak 
magnetic field but for all values of the interaction 

B. Quantum XY model 

In this section, we show that the attractive Hubbard model in the strong-coupling limit reduces to the quantum 
XY model (except in the low-density limit). The effective action of the phase 8r = —qr ~ —^r -I- Q • r of the 
superconducting order parameter is obtained by integrating out the p field. To quadratic order in p, the action reads 



s[p,e] 



-po'dr — i sin OoPrOr + 2hQ cos 9op^ 



J 



C0s(26'o) ^ {Pr + Pr 



(r,r') 



[l - cosiQr - er')][sin^ Bo + sm{2eo){pr + Pr') + 4:Cos^ OoPrPr' - 2sm^ Ooipl + pD] y (97) 



(r,r') 



Since the fluctuations of p are small, we can neglect all (i.e. Oq ^ n), which corresponds to the low-density limit 



the terms but sin in the coeflicient of 1 — cos(9r 
6r'). This approximation breaks down when sin 6*0 



of the attractive Hubbard model. Integrating out the p 
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field, we obtain the action of the quantum XY model 



s[e] 



2^°^^ + 167 



J2 [i-cos(er-er0]|.(98) 

(r,r'> ^ 



Taking the continuum limit, we recover the phase-only 
action derived in Sec.|ll|for C/ > 4t [Eqs. H52l58|l ]. 



C. Low-density limit: Gross-Pitaevskii equation 

In the strong-coupling limit, fermions form tightly 
bound pairs which behave as hard-core bosons (the hard- 
core constraint comes from the Pauli principle which pre- 
vents double occupancy of a lattice site). In the low- 
density limit, the hard-core constraint does not matter 
anymore and we expect to recover the usual action of a 
Bose superfluid, and in turn the Gross-Pitaevskii equa- 
tion. 

Using Eqs. II47I48II . 



(-lr[pr(2-pr)]'/'eT^''^ 



f^r = Pr-1, (99) 

we deduce from the Heisenberg model [Eq. Ij89|l ] 



S[p, 6] = J dT j d?r\^ -prOr + Jpl - {ho + 2J)pr 



J 

+ 4 



4pr 



Pr(Ver 



(100) 



in the low-density limit {pr <C 1) and to second order 
in gradient. We have taken the continuum limit in real 
space. Introducing the boson density pbr — Pr/2, we 
recover the action of a Bose superfiuid, 



(Vphr)^ term and replace pbr by pbo in the (V0r)^ term, 
we recover the action S[p, 6] (for x ^ 1) derived in 

secnini 

The classical equations of motion derived from the ac- 
tion (llfllll yield the Gross-Pitaevskii equation 



dr - Mfc 



2mi, 



*r+5|*r|'*r-0 



for the complex order parameter = 



(103) 



As 



pointed out in Sec. Ill C 21 equals the superconduct- 
ing order parameter Ar = | Ar|e*®'' in the strong-coupling 
low-density limit of the attractive Hubbard model. Alter- 
natively, the Gross-Pitaevskii equation can be obtained 
from the semiclassical spin dynamics. 



J 



r+S X Jlr — 2ho X fl^ 



(104) 



within a second-order gradient expansion and in the low- 
density Hmit Pr = f^r + 1 1- The sum over 6 in H104|l 
denotes a sum over nearest neighbors. 

For continuum models, the time- independent Gross- 
Pitaevskii equation has been obtained as the strong- 
coupling limit of the Bogoliubov-de Gennes equations 
It has also been shown that the results obtained in the 
strong-coupHng Hmit of an RPA calculation about the 
BCS state can be reproduced from the linearized ver- 
sion of the time-dependent Gross-Pitaevskii equationi^ 
In this section, we have directly shown the equivalence, 
in the low-density limit of the lattice case, between the 
low-energy e ffectiv e action of the superfiuid order param- 
eter Ar = ■y/pr/Se*®"' and the action of a Bose super- 
fiuid. Our approach is not limited to the low-density 
limit. Eq. 1110411 . together with (l99ll . holds for any density 
and can be considered as a generalization of the Gross- 
Pitaevskii equation to the lattice case. 



s[pb, e] 



where ^'r 



dr I d"^r^ ipbr^r + |/Ofcr - ^^bPbr 



1 



2rab 
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dr 



4p6r 







+|(«':*r)' 



1 



nib = 



P6r(Ver)' 



( 9r - p - ^ I *r 



2nib 



J' 

P6 = 2/10 + 4J, 
9 = 8J, 



(101) 



(102) 



are the mass, chemical potential and interaction constant 
of the bosons, respectively. Note that if we neglect the 



V. CONCLUDING REMARKS 

In this paper, we have studied the 2D attractive Hub- 
bard model using the mapping onto the half-filled repul- 
sive model in a uniform magnetic field coupled to the 
fermion spins. Our approach reproduces, in a unique 
framework, a number of previously known results. 

One of our main new results is the derivation of a 
low-energy effective action S[p, 8] which is valid for all 
values of doping a: = 1 — po and interaction strength 
U. S[p, Q] has been obtained by integrating out ampli- 
tude fiuctuations (|A|) and therefore does not describe 
the full dynamics of the superconducting order parameter 
Ar — \Ar\e^^'. Nevertheless, it is similar to the action of 
a Bose superfiuid with order parameter ^r — \J Pr/2e'®"' 
where the mass of the "bosons" and their mutual inter- 
action depend on x and V . This ensures that a Fermi 
superfiuid, as described by the 2D attractive Hubbard 
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model, will behave similarly to a Bose superfluid and ex- 
hibit the same macroscopic quantum phenomena. The ef- 
fective action S[p, 6] also describes the smooth crossover 
between the weak-coupling BCS limit and the Bose limit 
of preformed (local) pairs. 

Another important result obtained by our approach is 
a complete description of the phase diagram of the 2D 
attractive Hubbard model. From the phase-only action 
S[Q], we are able to extract an effective XY model and 
in turn the BKT phase transition temperature Tbkt- 
We identify a RC regime of superconducting fluctua- 
tions in the temperature range Tbkt <T< Tx and an 
incoherent-pair regime (with no superconducting short- 
range order) for Tx < T < Tpaii - The values obtained 
for Tpair and Tx are in good agreement with numerical 
results (when available). Near half-filling, we find that 
Tbkt is suppressed due to the strong q = (tt, tt) charge 
fluctuations which enlarge the symmetry of the order pa- 



rameter to SO (3). 

In the strong-coupling limit, the attractive Hubbard 
model maps onto the Heisenberg model in a uniform fleld. 
The latter reduces to the quantum XY model (except 
for a weak field, i.e. in the low-density limit of the at- 
tractive model). In the low-density limit, we recover the 
usual action of a Bose superfluid (including the terms 
proportional to (Vpr)^) and in turn the Gross-Pitaevskii 
equation. 
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APPENDIX A: HF CURRENT-CURRENT CORRELATION FUNCTION 



In this appendix we calculate the HF current-current correlation function 

K''f^'{q,q')^{j;{q)j;'i~'i'))c, (Ai) 

for q,q' — 0, Q. The calculation is performed at finite temperature in the classical limit: lim^^o = liniq^o linit^^^o 
and limqj^Q = limq^q limt^^^o- is the Fourier transformed field of 

f^r = -lt^l<j''(j)r+(, +C.C. {fi^x,y). (A2) 



We have 

JM(?~) = 7^E-M(k,q)0i-''0,+,, 

k 

vo{k,q) = 1, 

v^^oiK q) = -ii(e'('=-+'-) - e-'*^-). (A3) 

k — (k, icu) where oj — nT{2n + 1) (n integer) is a fermionic Matsubara frequency. In |IA1|I . one must have q — q' or 
q = q' + Q. The correlation functions of interest are 



n;:;:'(0,0) = ^ [z;,.(k,0)v(k,0)<.,<.,G.,(fc)G.,(fc) 

+v^,{k, 0) V (k + Q, 0)a^^,,a^;^^F,, (fc)F,, (fc)] , 

n;;;::(Q,Q) = -75^ E K(k,Q)v(k + Q,-Q)<.,<.,G.,(fc)G.,(fc + Q) 



k,(Ji,cr2 



+v^(k, Q) v(k, -QX,,,a^;^,T,, (fc)F,, (fc + Q)] , 
K^^^'iOM) = E K(k,0)v(k,-QX.,a^:^,T.,(fc)G.,(^) 



+i;^(k,0)v(k + Q,-Q)<,,<.,G.,(fc)T,,(fc)], (A4) 
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where 



w^(k, 0) = 5^,0 + (^^#02* sin fc^, 

t;^(k + Q,0) = 7M«M(k,0), 

u^(k, Q) = Sf^,o + 6f_,^o2itcoskf^, 

w,,(k + Q,-Q) = 7^v^(k,Q), 



(A5) 



and 7p = 5^,o — S^^q. We use the notation Sf^^o = 1 ^ (5^,0 and cr — —a. The HF propagators G and F are defined 
in ifTfiji . We have 



^^G,(k,ic^)G,,(k' + Q,^ 

LJ 

i^F<,(k,zc^)i^,,(k',^ 
i^F,(k,ic.)G,,(k',. 





£^ka H 


ekCTCk'tr' 
-Eko- 


Tk. ( 




ekfjek'CT' 

Sk. 


/Tk, 


Tk' 




V^^k^ 


Sk 




/Tka 


Tk' 






Fk 





Tk'a' £;k^ 



ekCTCk'tr' 
£^k'cr' 

ekCTCk'CT' 

i^k's' 



(A6) 



where Tkcr = tanh(/3£'kcr/2). Performing the sum over ci, tT2 in ljA4|l . we find that the only non- vanishing correlation 
functions are [using H);;);', (q, g') = Tiy^{q',q)\ 



nr,(o,o) 
n^^(o,o) 
n;i^^(o,o) 



no;;^o(o,Q) 



«'(k,0) 

<(k,0) 
Eh - El 



7k| . Hc;2 

kT 

TkT 



k -C/kT - -c^ki L-C'kT 



«2(k,0) 



k -*^kT ^ -^ki 



2^ %^z;.2 
ki 

2^ ^ki,„2 



^^{El, + .kT^ki - 7mA;?^^) - ^^{El, + .kT.ki - 7.A«s^) 



i;,.(k,Q)v(k,Q) 

-^kT -^ki 



'kl 



nZ,(Q,Q) 



^ ( -2 + ^(1 + 7.) 1 " C/kT^(l + 7.) 



T, 



Fi 



ill 

'kT 



.k l^kT^ + C/kt^ 



nSo(o,Q) = Af l^ekT (||-|^ 



ekCkT 



Tkt ?7kT 
, -^kT ^kT . 



n5„^(0,Q) = 2Af 



/TkT Tki 

k -E'kt - -E'kx V-E^kT -Bkx 



AHS 



AHS2 



-2+-^(1-7m) I -C/kT^(l-7M) 



^kT 



(A7) 
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where Uka — /3/2 cosh^ ^^^^j and ^p^^' — (5^,oi5^i',o + Sf^^o^i-i'^o- In order to obtain Eqs. (IA7|I . we have used the 

gap equations (I22I23|I and the symmetry relations eko- = — ek+Qff, -^ko- = -Eic+Qct, w^5^o(k, 0) — — Wp^o(^k, 0) — 
-i'^#o(k + Q,0), etc. 



1. Strong-coupling limit (/ > 4f (T = 0) 
Expanding Eqs. ijATII to leading order in l/U, we obtain (T = 0) 

n^o(0,0) = |sin2 0o, 

n§o"(Q,Q) = Icos^^o, 



oj-2 

nS^(Q,Q) - _cos0o(2-3cos2^o), 
noo(0,Q) - -^sin(20o), 

OJ.2 

no^(0,Q) = -^sinf?o(l-3cos2 0o). (A8) 

Here we consider only the correlation functions that are useful for the derivation of the effective action S[p, A] and 
the calculation of the velocity c of the phase collective mode. 



2. Correlation functions for ho = and fl^' = ( (T — 0) 

In Sec. mil we need the HF current-current correlation function at half filling and for a magnetization parallel to 
the z axis. They can be deduced from Eqs. ljA7ll with ho = h = and by making the rotation in spin space x ^ z, 
z — > y and y ^ x. The only non- vanishing correlation functions are then (the notations are those of Sec. Illlll 



nSo"(o,o) 




nr. (0,0) 




n^^(Q,Q) 


= wy^{QM) = - j^^J.S,.oEl 


nr^(Q,Q) 




nr.^o(o,Q) 


mo f el 


no,',^o(o,Q) 


mo f el 
2 LeI' 



(A9) 



where = a/c^ + 



APPENDIX B: EFFECTIVE ACTION S[m,n] 



In this appendix, we derive the effective action S[m, fl] = S[p, q, m] [Eq. lIMjl ] of spin fluctuations in the repulsive 
Hubbard model. 
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It is convenient to express the rotation matrix i?r as Rr = M(J7r)Af^(r2^') where M(ftr) is defined by 
M{nr)(T^M'<{nr) = cr • r2r: 



We then obtain 



^Sr = (-1) 



cos(f)e-2¥" _sin(^) 



sin(%) 63'' 



cos e 2 



9r + I sin(2pr) — sin 0oZ(5to 



HS 



— 



- ( + ^0 ) cos(2pr) + ^0 " cosOoiSmf^, 



A'i = COS 



Qr - Qr' 



Al^, = -z(-l)'-sin 



Al. = ^{-ly 



At „/ 



cos (pr + Pr' ) ^ 1 
2 



sin (pr - Pr') : 
sin (pr +Pr') ) 



cos (pr -Pr') • 



(Bl) 



(B2) 



Eqs. IIB2|I can be rewritten as 



Air = (-1)' 



-Tin 



4° 



2 (Pr<7r - Pr9r) + 2/ioPr " sin^oifJmP 



—q^ + 2ft,oPr — cos OoiSm 



HS 
r ' 



2 8 

:(-!)'' (gr - 9r')(Pr " Pr' ) > 



(B3) 



to quadratic order in pr, (5m^^, pr, gr and — 9r'- 

Keeping terms up to second-order in pr, (5m^^, Srrir, Pr, 9r and gr — qr' (r, r' nearest neighbors), the effective action 
is given by first- and second-order cumulants of Si and S2 with respect to the HF action: 



1 r'^ 

S[p,q,w^^M = {Si + S2) --{{Si+S2f)c + J drY, 



5ml + i5m^'^5mr - ^m^^Sm^^ 



2 



(B4) 



Using 



- 1, 



ia-c^r) = (-l)--2Ao, 
iaVr) - 0, 



(B5) 



where the averages are taken with the HF action, we obtain 

rl3 



{Si 



drY: 



iAo(pr(7r — Prqr) + 4Ao/ioPr ~ 2Ao sin 6*0 j(5m^^ + 2^qr ^ "^xhopl + X cos 9QiSm^^ 



(B6) 
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Similarly, from (r, r' are nearest neighbors) 



At 



>"'/'r') = 0, 



= 0, 



we deduce 



1-13 



(Pr + Pr 



2t ' 

x2 , ('J'r - ^r' 



Jo 



(B7) 



(B8) 



We have introduced the mean value {K) of the kinetic energy per site in the HF state. 

When calculating the second-order cumulant, it is sufficient to consider ^or and Ar,r' to linear order in Pr, 5m^^ , 

Pr, Qr, and Qr - Qr' : 



AX 

^Or 


= {-iy{2hoPr 

= -^i-iypr, 


— sin 


A^ 
^Or 


i 

— Or — COS 

2 


9oiSm 


^r.r' 

AX 

^r.r' 
Air' 


- 0, 
= 0, 

= i{-inpr + 


Pr'), 









HS 



In order to evaluate ({Si + 82)^)0, we write the action 5*1 + 5*2 as 



S1 + S2 ^ / dT V A] 
Jo 



ly(tot) ■I' 
fir Jfir^ 



= , X , y 



where A'j^j. — ^r.r+/i fo^' A* = 2;, The spin-density current j^^ is defined by l|A2ll . We then have 



S2)\^ f drdr' <*°'Hr)n-;(r,r;r',r')<,l.^°*)(r'), 
"'0 



where 



U;';,,{r,T;r',r') = U;^{T)j;,,,{T')), 



is the HF current-current correlation function. IIJ^^, is calculated in Appendix El 

In order to calculate the second-order cumulant, we write the actions Si and S2 [Eqs. ll32ll ] as 



OrJOrJ 



51 = f dr ^Or 

52 = -i dT Y ^^rJ^r. 

Jo 



(B9) 

(BIO) 
(Bll) 

(B12) 



(B13) 



(B14) 
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where A^^. = M = Aq^. and AJ;^^. are given by Eqs. (IB9|I . The spin-density current j^^ is defined in 

ljA2|l . In Fourier space, we obtain 



where 

We deduce 

((^1+^2)') 



^0(9) = -(^i^, 



(B15) 



(B16) 



E E 

1 



i-q - Q)z^. (g' + Q)p-qPq'Ull, {q + Q,q' + Q) + -z^{-q)z^, {q')q-qqq,W;^, (9, q') 

1 



- E 

fi—0,x,y 



zM' + ^V-qPq'^Z^'l + Q, 9' + Q) + -z^{q')r^qqq,I{ll(q + Q, q') 



+z 



i-q - Q)p-qrq,nlUq + Q, q + Q) + ^^^(-^I'z-g-'^g'n^Slg, q + Q) 

z^.i-q - Q)p-qiSmfU%iq + Q, q') + ^z,,i-q)q^qiSmfli;i{q, q' 
r-qrq,nZ^,{q + Q, g' + Q) - cos^ e^5vnl\6mfW,^^{q, q') 

COS 00 [r-qi5mflil^{q + Q, g') + z5mH?r,.ng5(q, g' + Q)] I , 



(B17) 



where = 2hoPq — sin OoiSrn^^ . To proceed further, we use the expression of the correlation function Il'j^''^,{q,q') 
obtained in Appendix El We thus have q — q' in ijBlTII . In order to obtain the effective action of the q field to 
lowest order in d^i,q, it is sufficient to retain the first-order derivative terms Since = 0{d^), one can use 

zo(g) — —uju, z^^o{q) = —iqfj., za{q -f Q) = 0, z^^o(9 + Q) = 2, and evaluate the correlation functions HJ;^^/ at 
q g' = in ijBlTII . This gives 



{{Si 



S2?) 



J2\p-9p4^ E n^^(Q,Q) + 4/jXo(Q,Q) + monS,^(Q,Q) 



+g_^g^(-^5in-(o,o) + ^n-(o,o)). 



-Sm'^pmf (sin2 0onSo"(Q,Q) +cos2 0ongg(O,O) -Hsin(20o)ngg(O,Q)) 
- [p-qqqco. (2n^^(0, Q) + /iongS(0, Q)) - CO.] 

-i[p_g5mf (2/10 sin 0onSo"(Q, Q) + 4sin(?ong^(Q, Q) 
+4 cos OoU'^l (0, Q) + 2ho cos 6»ongg (0, Q)) + c.c. 



q-gSmf"^ (sin0ongo"(O,Q) +cos0ongo(O,O)) - c.c. 
A comment is in order here. The first-order cumulant IIBfi|l gives a term linear in 

6Sb = ^a; / dr^q^. 



(B18) 



(B19) 
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This term comes from the Berry phase term Sb = dT^^{(l)\.R\Rr4)r) ■ The SU{2) matrix i?r is defined up to 

the U(l) gauge transformation Rre~i^''^'^i\ The Berry phase term depends on the gauge choicei^ The 

gauge-dependent term is given by 



1^ dTy^{4>la-nfd,Mr = -l:nio drY^i^,. (B20) 

In the following we take -i/^r = 'Pr/'mo- As shown in Sec. Ill C\ this choice ensures that in the attractive model, half 
the fermion density is identified as the conjugate variable of the phase of the superconducting order parameter, as 
required by gauge invariancei^ Including the gauge-dependent term of Sb [Eq. IIB2fl|l ] in 5Sb, we obtain 

/•/3 

5Sb = -1:Po / dr^gr- (B21) 







From Eqs. l|B6IB8IB18IB21ll . we then obtain 

+2p^qUp^Hs{q)6mf^ + 2q^qUg„,Hs{q)Smf^ + 2iSm,^qSmf^ - ^Sm^qSmq] + 5Sb, (B22) 

where 

npp(9) = -4 E n^"^,(Q,Q)-4/i2ngo-(Q,Q)-16/ionS^(Q,Q)-4a;/io + 4(-if), 

u,,iq) ^ ^n-(o,o) + ^(-i^), 

n„Hs„Hs(g) = sin2 0ongo-(Q,Q) + cos2eong;i(O,O)+sin(20o)ngo(O,Q), 
np,(g) = c^,[2n^^(0,Q) + /iongg(0,Q) + 2Ao], 
n^^Hs (g) = 2i [ho sin ^ongo" (Q, Q) + 2 sin eoU^'^iQ, Q) + 2 cos doU'ol (0, Q) + /^o cos 0ongJ(O, Q)] , 

n,™Hs(g) = ^c^,[sin0ongg(O,Q) + cos(?ongo(O,O)]. (B23) 

Pg, (7g, (5m?^ and (5™^ are the Fourier transformed fields of Pr, 9r, Sm^^ and (5mr. = {q,iuj^) where ojiy = v2t:T 
{u integer) is a bosonic Matsubara frequency. Eqs. llB23p are valid in the hydrodynamic regime (g ^ 0). They are 
sufficient to obtain the effective action of the AF field q to order 0{df^). 

Integrating out the Hubbard-Stratonovich field m^^, we obtain Eq. (j34|l with 

flgmHS n^HSg 

J-7-V) HS jyi HS 



n,™ - (B24) 

The action S[p, q, m] takes a very simple form in the weak-coupling (Slater) and strong-coupHng (Mott-Heisenberg) 
limits. 
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Slater limit 



We assume that we are not too close to half-filling so that the zero-temperature order parameter A[/^ is given by 
26)l . Since Aq^ is exponentially small at weak coupHng, Oq ~ n and irn^^ ~ xU/2 [see Eqs. II15I21II ]. We also have 



n§S(o,Q) 



aHS 
^0 



4t 



de- 



ep 



° '-4t [(e-ei.)2 + AHs2]3/2 

~ 0. (B25) 

Since the integral in ljB25ll is peaked around e — ep for A^^ 0, we have replaced the density of states A/'o(e) by its 
value at the Fermi energy ep [see the discussion after Eq. Ij25|l ]. For the same reason, we can extend the integration 
range to ] — oo,oo[, so that the integral vanishes. This result is a consequence of the particle-hole symmetry which 
holds in the weak-coupHng (BCS) limit of the attractive model. Similarly, we find ng^(0, Q) ~ 0. We therefore have 



n^^HS^HS 



Ilg^HS 



n§g(o,o) 



ngo(o,o), 

2w^Ao ~ 0, 



(?) = 0, 
(9) - - 



n-(o,oK. 



(B26) 



Using Eqs. ljB24ll . we deduce 



1 



ngg(o,o) 

~2 



U 
2"' 



(B27) 



and Ilpq — Ilpm — 0. p fiuctuations do not couple to q and m fiuctuations. We therefore obtain the effective action 



S[q,i 



dT 



d^r 



(-K) 



1 



u 



n-(o,o) 2 



^ 5 



8 

- SSb, 



(B28) 



where we have taken the continuum limit in real space. 



2. Mott-Heisenberg limit 

In the strong-coupling limit, there are well-defined local moments with a fixed amplitude so that we can ignore 
the fluctuations of m^^ and to. Low-energy fluctuations correspond to direction fluctuations of these local moments. 
This can be seen explicitly by integrating out the to fleld in the action S[p, q,m^^ ,m] [Eq. ljB22ll ]. This yields the 
replacement 

/ 2 

n^HSmHS n^HS^HS n„jHSmHS - — . (B29) 

To leading order in 1/U , we have [see Eqs. ljA8|l in Appendix fX] 

2 

n^Hs„Hs(g) = - — , 

IIpmHS = O I — 



C/2 

IP 



n,™Hs = o ( ^ ) . (B30) 
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Integrating out yields the correction terms 



IT 

nqm'^^ n^HS q 

n' 

^HS^HS 

n^m^^ n^Hs ^ 



= o 



o 



C/3 



o 



C/3 



to Tlpp, Tlqq and Ilpg, respectively. These terms can be ignored in the limit U ^ 4t. 
Using Eqs. IIA8|I of Annendix IA II we therefore have 

^ppi<j) = 8Jsin^6'o, 

^qq{q) = q^^sin^fi'o, 
npg(g) = u;,. sin 6*0, 

and Upm^s = Hqm^s ~ 0. Using then II = 11 [Eqs. (IB24|l ]. we obtain the effective action 



S[p,q] = - dr d\ 



1i\J\ — x'^PrQr 



-j{l-x^){Vq.r + 8J{l-x')pl 



SSf 



APPENDIX C: KINETIC ENERGY IN THE HF STATE {T = 0) 

In this Appendix, we derive Eq. ll56|l . We start from 

- Vk • Vk + Vk 



Vk.(flvk) ^ fIVk.vk + Vk(fl). 



Vk 



(B31) 



(B32) 



(B33) 



(CI) 



where Vk = Vk^k- Integrating the left hand side of this equation over the entire Brillouin zone, we obtain 

-Vk • dlk = 0, (C2) 



where the contour in the last integral is given by the Brillouin zone (BZ) boundary and dlk is perpendicular to the 
contour. The integral vanishes since Vk • dlk = at the BZ boundary. We deduce from l|(ni(::2|l 



SkT 



k £'kT 



(C3) 



In the weak-coupling limit, using Ag^^/E'^^ = 2(5(ek — ef), we obtain lf5fill . 



APPENDIX D: CHARGE {pr) AND PAIRING (A^) FIELDS 

In the section, we relate the charge (pr) and pairing (Ar) fluctuations of the attractive model to the fields p, q and 
5m defined in the repulsive model. We rewrite Sj [Eq. l(35ll ] as 



r0 

■Jo 



(Dl) 
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where Bor = i?J Jr ■ (rRr = Y.u=x,y,z Bor'^" and 

Bor = (-l)"" COs(0r-0o) [cos </JrJr"+ sin (/JrJrl -(-!)"" Sin(0r-0o)^r': 

Bl = -(-l)""sin^rJr +(-!)"" cos</PrJ;', 

= sin(0r - Oo) [cos (^r + siii if^jy] + cos(6'r " 6*0)7;. (D2) 

To first order in J and second-order inp, 5m^^,Sm, the effective action is given by {Si+S2+Sj)-'^{{Si+S2+Sj)'^)c. 
The first-order cumulant gives the source-dependent contribution 



Jo 



(D3) 



From the second-order cumulant, we obtain to linear order in the source 
Sf = -fdrdr' Bo^r(r)n„7;(r,r;r',r')A;:(r)(r') 

fj,' =Q,x,y 

= - [ «'^E{(-l)"'^Or [4nS^(Q, Q)Pr + nSo-(Q, Q)(2Mr - sin0oi<5mf ) + ng^(0, Q)(-^'7r - cos^oi'JmHS 

[AlilliO, Q)p, + ngg(0, Q)(2/ioPr - sin^o^^mP) + ng^(0, 0) (-^gr - cos 0oi<5mHs)] |, (D4) 

where the last expression is valid in the hydrodynamic regime. A'^'"^°^^ is defined in ijBllll . In the presence of the 
source J, the effective action S[p, q, m^^, m] [Eq. IIB22ll ] should be supplemented with S^j'' + S'^P . The integration of 
the m^^ field then leads to the source-dependent action + + s'f^ where — S^lsm^s^Q and 

Sf = -— ^ rdTVf(-iri3o-r[sinWo(Q,Q)+cos^ongS(0,Q)] 

+So\ [sin doKo (0, Q) + cos 60^0 (0, 0)] } 

We have introduced H^hs^ — ^musq/uji,. In Eq. (|D5ll . the correlation functions IlmHSmHs ,nmHSp,n^^^Hs^ are evaluated 
at 9 = 0. Taking the functional derivative of s'^P + + S^f"^ with respect to J, we finally obtain 

6pr = Prh4Ao-2/iongg(0,Q)-4n^^(0,Q)]+n-(0,0)^gr 

sin6longg(0, Q) -I- cos6'ongg(0,0)] (H^nSpPr + H^^Hs^igr + iSm^), 



n^HSy^HS 



5|Ar| = pr[-x-/iongJ(Q,Q)-2ng,^(Q,Q)]-f ng5(0,Q)^gr 



[sin0OnSJ(Q, Q) + COS0ongo(O, Q)] i^m^^pPr + ^In^Sqiqr + iSm,), 



Or = -9r. (D6) 

1. Strong-coupling limit 

In the strong-coupling limit, we can obtain a simple relation between pr, and rrirflr- In Sec. II VL we show that 
the Hubbard model reduces to the Heisenberg model when U ^ At. We can carry out the same derivation in the 
presence of the source term ljDl|l . Integrating out the fermions, we obtain the source-dependent term 

Sj= / dTy^{cblRlJr-<TR,(br)e.t (D7) 
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to leading order in 1/U . The average in ljD7|l is taken with respect to the atomic action H8(i|l . The matrix i?r satisfies 
Rrd^Rl = o- • rir (see Sec. lIV Ajl . Using (jir) at = Su,z and i?JJr ■ crRr = cr ■ TZ^^Jr, where 

cos 6r cos (pr — sin (fir sin 9r cos (fir 

TZr = I cos 9r sin </?r COS </?r sin 9r sin (^r I (D8) 
— sin 9r COS ( 

is the so (3) rotation matrix which maps z onto fir, we obtain 

r/3 



Sj — dr Jr • fir 



(D9) 



From Eqs. II36|I we then deduce 



2 



where = n^±iny. 



APPENDIX E: HF CUMULANTS {Sp + + Ssm) AND {{Sho + Ss,r^f + 2iSh„ + 5'i„)(5p + Si))c 

In this appendix, we calculate {Sp + Sho + Ssm) and {{Sho + Ssm)'^ + 2{Sho + Ss,n){Sp + Si))c (Sec. IIIip . The 
first-order cumulants read 

rl3 



Jo 

i— ^ / dr V'(-l)'"(.^rC0s6'r + V'r), 



Jo 

[' drY^i-DVr, (El) 
Jo 



u 



and 

{Ssn 



drJ2i-iySmr{j^r) 

Jo 

/ rfrVtJmr, (E2) 

"'0 



2too 



where we have used — cos^r and the saddle-point equation 2mo/U = (— l)''(c]!(7^Cr). Since iir is slowly varying, 
{Sho) vanishes. From HE2I69|I . we conclude that there is no Hnear contribution in Sm. 
Let us now consider the second-order cumulant 

(Sh,Ssm) = ho / drdr' V B^^iT)Il'^^ir,T;v' ,T'){-iy' Smr> 

= ho B-oi~mooiq,'i')S^g'+Q- (E3) 



1 _ 

(— l)''z). It is given in Annendix lA 21 Since Bq^. is slowly varying, we can evaluate Il'^'^,{q,q') at g = in order to 



nj^^, is the HF current-current correlation function for hg — and an AF order parallel to the z axis (i.e. $7' 
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obtain the result to second order in ho, and 5m. Since ^Q^iq = 0,q') = f Appendix lA 2|l . (ShoSsm) vanishes. A 
similar calculation shows that {SpSsm) = {SiSsm) = 0. We therefore conclude that amplitude fluctuations decouple 
in the limit of a weak magnetic fleld. 
The contribution due to S*^^ is given by 

(^^o) - hi f drdr' Bo"rMn[;o^'(r,T;r',T')<,(r') 

u^u' = x,y,z 

= K B-{-q)Woo{q,q')BI^'{q'). (E4) 



Since Bq^, is slowly varying, we can evaluate Hon at g = = 0: 

(52J=/jXo Tdr^sin^^r. (E5) 

Here and in the following, we use the notation Uqq = Uqq {q — 0,q' — 0). We have used IIoq cx Si,^^i{6^^x + ^v,y)- A 
similar calculation gives 

{ShoSp) = /longiJ / dr V-sin2 6'r93r, 

Jo r 

(ShoSi) = homoU^oS rfrVsin^rl-cosV-r/^ + sinV-r^i^). (E6) 
Jo 

To express (5^^), {ShaSp) and {ShgSi) in terms of n and 1, we use 

(ho X rir)^ = hi sin^ 9^, 
ho • (rir X hr) = hlipr sin^ 9r, 

= sin6'r(- cos^Ar'? + sin-^r^r)- (E7) 

The last result follows from Lr = TZAr and 1^—0, where TZr is given by 

cos dj. cos (fr COS Ipr — sin Ifir sin ^/^r — cos 9r COS Ifi^ sin Ipr — sin 9?r cos Ipr sin cos (fr 
TZy. — I COS Or sin cos '(/'r + COS ipr Sin V'r — cos Or Sin sin V'r + cos (/?r cos f/jr sin 6'r sin ifr I . (E8) 

— sin Or COS ■i/'r sin Or sin ^/^r cos 



We therefore obtain 



Jo 

(ShA) = ^ngo^l dr^ho-Kxnr), 



= moHgo" / drVho-Lr. (E9) 
Jo 
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